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FOREWORD 


This  report  completes  a  broad  research  pro¬ 
gram  in  newconcepts  for  guidance  and  control  of  tactical 
missiles .  Other  work  performed  under  this  contract  is 
described  in:  TR-17G-1,  ’'Adaptive  Control  and  Guidance 
for  Tactical  Missiles”;  TR-170-2,  "Optimal  Stochastic 
Guidance  Laws  for  Tactical  Missiles”;  TR-170-3, 
"Adaptive  Control  with  Explicit  Parameter  Identifica¬ 
tion  for  Tactical  Missiles”.  The  authors  wish  to  express 
appreciation  to  Mr.  David  Siegel  of  the  Office  of  Naval 
Research  for  his  encouragement  and  support  throughout 
this  investigation. 
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1. 


INTRODUCTION 


1.1  BACKGROUND 

Previous  investigations  of  homing  guidance  laws  for  tactical 
missiles  have  discussed  a  number  of  guidance  techniques  derived  using 
principles  of  optimal  stochastic  control  theory  (Refs.  1-5), which  offer 
potential  improvement  over  conventional  proportional  guidance.  In  the  past, 
comparative  evaluations  of  these  guidance  laws  obtained  from  computer 
simulations  have  generally  been  based  upon  simplified  linear  models  of 
guidance  system  dynamics.  It  is  found  that  certain  frequently  neglected 
effects,  such  as  limited  missile  maneuvering  capability  and  dynamic 
coupling  between  the  missile  airframe  and  its  homing  seeker,  signifi¬ 
cantly  influence  missile  performance.  Although  the  simplified  models 
are  adequate  for  obtaining  qualitative  comparisons  of  performance  for 
different  guidance  laws,  definitive  quantitative  evaluations  require  more 
accurate  mathematical  representations  of  the  missile -target  engage¬ 
ment  situation. 

The  purpose  of  this  study  is  to  provide  comparative  perfor¬ 
mance  evaluations  of  several  modern  guidance  laws  using  a  computer 
simulation  based  upon  a  mathematical  "truth"  model  of  the  guidance 
system.  The  model  includes  descriptions  of  the  following  critical  and 
often  neglected  missile  characteristics:  autopilot  dynamics,  acceleration 
limit,  homing  seeker  dynamics,  range  dependent  seeker  noise,  geomet¬ 
ric  nonlinearities,  and  random  target  motion.  The  guidance  laws  to  be 
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compared  include  proportional  guidance,  optimal  linear  guidance  and  op 
timal  nonlinear  guidance.  These  laws  are  derived  based  upon  a  simpli¬ 
fied  "design"  model  that  neglects  many  of  the  above  effects;  then  their 
performance  is  tested  in  a  simulation  of  the  truth  model.  This  deter¬ 
mines  the  sensitivity  of  the  guidance  system  to  errors  in  the  design 
assumptions. 


1.2  METHODOLCXJY 

Traditionally  the  analysis  of  random  effects  in  a  nonlinear 
system  is  accomplished  using  the  monte  carlo  simulation  technique. 

That  is,  a  large  number  of  sample  solutions  to  the  randomly  excited 
differential  equations  of  motion  are  generated  by  computer,  and  ensemble 
averages  of  the  variables  of  interest  are  computed.  Frequently  a  large 
number  of  trials  (e.  g. ,  several  hundred)  are  required  to  obtain  a  suffi¬ 
ciently  accurate  estimate  of  average  performance.  This  procedure  is 
quite  expensive  in  terms  of  computer  running  time  and,  for  a  complex 
system,  only  a  few  different  system  configurations  can  be  evaluated 
within  a  limited  computational  budget. 

Recently  a  new  approximate  analytical  technique  has  been 
developed  for  calculating  the  statistics  of  nonlinear  stochastic  systems 
(Ref.  6).  This  method,  called  the  Covariance  Analysis  Describing 
Function  Technique  (CADET),  yields  good  agreement  with  monte  carlo 
results  and  requires  significantly  less  computer  running  time.  The  basic 
principal  of  CADET  is  that  the  system  nonlinearities  are  linearized  using 
describing  function  theory;  then  the  statistics  of  the  resulting  linearized 
system  model  are  analyzed  using  covariance  analysis  techniques.  The 
saving  in  computer  time  arises  because  only  one  solution  to  a  matrix 
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differential  (or  difference)  equation  is  required  to  determine  the  statistics 
of  a  given  system  configuration,  as  compared  with  several  hundred  solu¬ 
tions  of  the  system  differential  (or  difference)  equations  with  the  monte 
carlo  method.  The  CADET  technique  is  attractive  for  analyzing  guidance 
systems  because  it  permits  the  investigation  of  several  different  guidance 
laws  over  a  wide  range  of  values  of  the  system  parameters.  Consequently, 
the  approach  taken  in  this  study  is  to  perform  most  of  the  guidance  sys¬ 
tem  performance  evaluations  with  a  CADET  computer  simulation,  validat¬ 
ing  the  latter  with  a  monte  carlo  simulation  in  a  selected  number  of 
cases. 

1.3  ORGANIZATION 

Chapter  2  of  the  report  outlines  the  principal  features  of  the 
simulation  truth  model  and  the  guidance  laws  under  comparison.  Chapter 
3  presents  simulation  results;  a  summary  and  list  of  the  major  conclu¬ 
sions  are  given  in  Chapter  4,  Appendix  A  provides  a  detailed  discussion 
of  the  system  truth  model;  the  theoretical  background  for  CADET  is  pre¬ 
sented  in  Appendix  B. 
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2.  GUIDANCE  SYSTEM  MODEL 

2.1  System  Dynamics 

The  principal  elements  of  a  missile  guidance  system  are  indi¬ 
cated  in  the  block  diagram  of  Fig.  2.1-1 .  This  section  briefly  describes 
the  function  of  each  subsystem,  exclusive  of  the  guidance  law,  and  the 
features  included  in  the  simulation  model.  The  missile  guidance  law  is 
treated  in  Section  2. 2.  The  reader  is  referred  to  Appendix  A  for  more 
complete  details. 


MISSILE 


DISTANCE 


Figure  2.1-1  Guidance  System  Truth  Model  with 
Major  Subsystems 
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Seeker-  The  homing  seeker  has  the  function  of  tracking  the 
target.  Fundamentally  it  measures  the  angular  position  of  the  line -of  - 
sight  (LOS)  to  the  target  relative  to  the  sensor  centerline  by  observing 
the  direction  of  electromagnetic  energy  reflected  from  the  target  (as  in 
a  radar  seeker)  or  transmitted  by  the  target  (in  the  case  of  an  infrared 
or  electrooptical  seeker).  This  measured  error  signal  generally  serves 
two  functions:  first,  in  the  case  of  a  gimballed  seeker,  the  error  is 
used  to  drive  the  seeker  antenna  so  as  to  maintain  a  small  tracking  error; 
this  is  necessary  to  prevent  losing  track  of  the  target.  Secondly,  the 
response  of  the  seeker— e.  g. ,  its  gimbal  angle  rate,  the  magnitude  of 
the  tracking  error,  etc,  —is  used  to  infer  the  LOS  angular  rate,  or  the 
LOS  angle  measured  with  respect  to  some  reference  orientation.  Either 
of  these  quantities  may  be  used  in  the  guidance  law  for  generating  mis¬ 
sile  acceleration  commands  which  direct  the  missile  toward  intercept. 

The  model  of  the  seeker  includes  the  antenna  pointing  dynamics, 
as  well  as  parasitic  coupling  with  the  missile  airframe.  The  latter  can 
arise  from  an  error  in  the  direction  of  the  LOS,  as  perceived  by  the 
seeker,  caused  by  aberration  of  the  electromagnetic  energy  as  it  passes 
through  the  protective  covering  (radome  in  the  case  of  a  radar)  of  the  mis¬ 
sile.  This  error  typically  depends  upon  the  missile  attitude;  therefore  the 
airframe  dynamics  are  coupled  to  the  seeker  measurement. 

Autopilot  -  The  autopilot  refers  to  the  missile  airframe 
dynamics,  together  with  its  stability  augmentation  system.  It  can  be 
viewed  as  a  servomechanism  designed  so  that  its  output  acceleration 
follows  the  input  acceleration  commands  as  closely  as  possible.  This 
study  uses  two  selectable  models  for  the  autopilot— one  is  a  first-order 
lag  which  provides  a  first  approximation  to  actual  airframe  response; 
the  other  is  a  ihird-order  transfer  function  which  more  closely  represents 
the  dynamics  of  an  actual  airframe.  In  addition,  the  autopilot  model 
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includes  one  of  the  most  significant  guidance  system  nonlinearities— an 
input  command  limiter  which  represents  the  limited  maneuvering  capa¬ 
bility  of  the  airframe. 

Missile -Target  Kinematics  -  The  model  of  the  missile-target 
kinematics  assumes  that  the  relative  closing  velocity  (range-rate)  is  con¬ 
stant  and  that  the  missile  and  target  accelerations  are  normal  to  their 
respective  velocity  vectors.  This  implies  that  both  velocities  have  con-  _ 
stant  magnitude  but  variable  direction.  The  components  of  missile  and 
target  acceleration  which  are  normal  to  the  line-of-sight,  and  which 
control  the  terminal  miss  distance,  are  represented  as  trigonometric 
functions  of  the  appropriate  orientation  angles.  This  constitutes  another 
important  nonlinearity  included  in  the  model. 

Measurement  Noise  -  The  noise  in  the  measurement  of  line- 
of-sight  angle  is  caused  by  a  number  of  different  error  sources,  some 
of  which  are  range  dependent.  Three  types  of  noise  are  included  in  the 
model:  effective  receiver  noise  (decreases  with  decreasing  range),  target 
angular  scintillation  noise  (increases  with  decreasing  range),  and  range 
independent  noise  caused  by  the  seeker  servo  system  and  possibly  by 
amplitude  fluctuations  in  the  received  signal. 

Target  Motion  -  The  guidance  system  model  includes  a  pro¬ 
vision  for  a  randomly  accelerating  target.  The  latter  is  representative 
of  maneuvers  that  are  constant  in  magnitude,  but  switch  sign  at  random 
times— sometimes  referred  to  as  "jinking”.  The  random  motion  is 
modeled  as  the  output  of  a  low  pass  filter  driven  by  white  noise:  the 
latter  has  the  same  autocorrelatibn  function  as  the  jinking  maneuver. 

Initial  Conditions  -  A  number  of  initial  conditions  are  specified 
in  the  system  model  which  determine  the  missile  trajectory  and  influence 
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terminal  accuracy.  These  include: 

•  Launch  Range 

e  Missile  Velocity 

•  Target  Velocity 

•  Missile  Heading  Error 

•  Closing  Velocity 

For  this  investigation,  most  of  the  above  quantities  were  assigned 
nominal  values  that  were  held  fixed  throughout  the  study. 

Mathematical  models  are  developed  in  Appendix  A  for  each 
of  the  subsystems  indicated  in  Fig.  2. 1-1.  Figure  2. 1-2  provides  a 
detailed  block  diagram  of  the  complete  guidance  system.  The  missile- 
target  motion  is  assumed  to  be  restricted  to  a  single  plane  and  the  effects 
of  gravity  and  aerodynamic  drag  are  neglected.  The  kinematic  variables— 

^m’  vm’  vt’  V  9a’  e’  ajj’  \  are  defined  in  Fig.  2.1-3.  The  quan¬ 
tity  kr  is  a  gain  representing  the  aberration  error  which  dynamically 
couples  the  missile  airframe  motion  to  the  homing  seeker.  The  transfer 
function  coefficients  rj,  At,  a^,  a2>  ag,  b^  bg,  andbg  are  assigned 
values  to  yield  realistic  dynamic  characteristics.  The  switches 
Sj  and  S2  are  positioned  according  to  the  type  of  guidance  law  being 
evaluated;  their  settings  are  explained  in  Section  2.2. 

One  of  the  more  important  autopilot  parameters  is  ag,  which 
combines  with  kr  in  determining  the  high  frequency  open  loop  gain  of  the 
parasitic  attitude  loop;  if  the  latter  is  excessively  large,  system  insta¬ 
bility  can  result.  The  value  of  ag  tends  to  increase  with  missile  altitude 
because  a  larger  missile  angle  of  attack,  and  hence  a  higher  transient 
pitch  rate,  is  needed  to  generate  a  given  acceleration  as  the  air  density 
decreases. 
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Figure  2.1-2  uses  Laplace  transform  notation  to  indicate 
input -output  relations  associated  with  linear  subsystems.  However,  the 
system  performance  cannot  be  analyzed  in  terms  of  this  notation  because 
of  the  nonlinearities,  a  time-varying  gain  (1/r^)  and  time-varying  rms 
seeker  measurement  noise.  Chapter  3  discusses  the  methodology  used 
to  determine  the  rms  terminal  miss  distance. 


2.2  MISSILE  GUIDANCE  LAWS 

The  motivation  for  this  study  is  partially  provided  by  previous 
work,  described  in  Ref.  7,  which  investigated  optimal  stochastic  guid¬ 
ance  laws  in  the  presence  of  missile  acceleration  limiting,  using  a 
limited  number  of  monte  carlo  simulations  and  a  simplified  system  model. 
It  was  determined  that  under  some  circumstances  an  optimal  nonlinear 
law  achieves  as  much  as  fifty  percent  reduction  in  the  miss  distance 
achieved  with  the  linear  law.  The  optimal  guidance  laws,  together  with 
classical  guidance  methods,  provide  a  hierarchy  of  techniques  that  are 
potentially  applicable  for  tactical  missiles,  and  which  need  evaluation 
using  a  more  realistic  system  model. 

Generally  speaking, the  guidance  law  is  thought  of  as  two  cas¬ 
caded  functions  —  filtering,  or  state  estimation,  and  control.  The  function 
of  the  filtering  operation  is  to  obtain  estimates  of  those  variables  needed 
to  mechanize  the  control  law.  The  latter  prescribes  the  acceleration 
command  according  to  a  policy  which  will  direct  the  missile  trajectory 
to  intercept  the  target.  The  methodology  used  for  designing  the  various 
laws  is  discussed  in  Section  A.  6 .  Five  different  laws  are  considered; 
their  distinguishing  characteristics  are  summarized  in  Table  2.2-1. 
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TABLE  2. 2-1 
MISSILE  GUIDANCE  LAWS 


Identifying 

Symbol 

Distinguishing  Characteristics 

Type  of 
Filter 

A 

Proportional  Guidance 

First-Order, 

Low-Pass 

B 

Proportional  Guidance 

Kalman 

C 

Optimal  Linear  Guidance; 

Accounts  for  Target  Maneuvers 

Kalman 

D 

Optimal  Linear  Guidance; 
Accoimts  for  Target  Maneuvers 
and  Missile  Autopilot  Dynamics 

Kalman 

E 

Optimal  Nonlinear  Guidance 
Accounts  for  Target  Maneuvers 
Missile  Autopilot  Dynamics, 
and  Missile  Airframe  Saturation 

Kalman 

Law  A  is  simply  conventional  proportional  guidance,  preceded 
by  a  low-pass  filter  to  suppress  measurement  noise  in  the  sensor  output, 

§ '.  For  this  case  switches  Sj  and  Sg  in  Fig.  2.1-2  are  both  in  position  1, 
and  the  guidance  law  is  mechanized  as  shown  in  Fig.  2.2-1,  where  n'  is 
a  specified  constant  navigation  ratio  and  v  is  the  missile -target  closing 

V 

velocity*  (range-rate). 

Guidance  law  B  uses  a  Kalman  filter  to  provide  an  optimal 
estimate  of  LOS  angular  rate,  together  with  the  same  proportional 


*It  is  assumed  that  the  errors  in  measuring  vc,  rt  ,  and  a.  are  negligible. 
This  is  a  reasonable  assumption  in  the  case  of  a  radar  homing  seeker. 

However,  range  and  range  rate  are  not  accurately  known  in  a  system 
using  an  infrared  seeker. 
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guidance  law  as  configuration  A.  A  detailed  diagram  of  the  filter  is 

shown  in  Fig.  A.  6-3;  it  estimates  the  three  state  variables  y^,  and  afc, 

assuming  that  both  range  (r.  )  and  missile  acceleration*  (a .)  are  known. 

tm  z  ^ 

These  estimates  of  y^  and  y^  are  combined  to  yield  an  estimate  of  fl 
according  to 


Vd 


(2.2-1) 


The  switches  Sj  and  S2  (Fig.  2.1-2)  are  both  in  position  2  and  Fig.  2.2-2 
illustrates  the  mechanization  of  the  guidance  law.  The  Kalman  filter  is 
inherently  a  digital  processor;  therefore  the  acceleration  command  is 
computed  discretely,  rather  than  continuously.  Observe  that  the 
estimate  of  target  acceleration  is  not  used  for  control  in  configura¬ 
tion  B. 

Guidance  Laws  C,  D  and  E  are  all  represented  by  the  diagram 
in  Fig.  2.2-3.  In  configurations  C  and  D,  the  control  laws  are  chosen  to 
minimize  the  performance  index 


=  lim  E  ]  (miss  distance)^  +  y  f 

V.*  n  \  J  ~ 


aj  (t)  dt 
c 


(2.2-2) 


making  use  of  known  results  from  optimal  stochastic  control  tneory. 
Configuration  C  minimizes  J,  neglecting  autopilot  dynamics,  resulting  in 
c^  =  0,  n '  -  constant,  and  Cg  is  determined  as  a  function  of  target  band¬ 
width  (A.)  and  time-to*  go  until  intercept  (t  ).  Configuration  D  is  derived 
including  a  first-order  model  of  autopilot  dynamics  having  bandwidth  Am; 
for  this  case  Cg  is  the  same  as  in  lawB,  whereas  n' and  c^  are  both  functions 


See  footnote  on  page  2-7 
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of  X  ana  t  .  Both  laws  C  and  D  are  optimal  linear  guidance  law&;  the 
m  go 

principal  difference  between  them  is  that  in  law  D,  n '  increases  with 
decreasing  range  and  becomes  very  large  near  intercept  to  compensate 
for  the  effective  lag  in  missile  acceleration  caused  by  the  dynamics  of  the 
missile -autopilot  combination. 


Configuration  E  is  designed  to  minimize 

J  =  E  1 1  miss  distance  1 1  (2.2-3) 

subject  to  the  explicit  constraint 


max 


(2.2-4) 


Because  the  limit  on  missile  acceleration  is  accounted  for  in  the  problem 
formulation,  the  resulting  optimal  nonlinear  guidance  law  differs  from 
those  described  previously.  The  mechanization  takes  the  same  form  as 
that  in  Fig.  2.2-3;  the  only  difference  is  the  manner  in  which  n'  is  com¬ 
puted.  Typically,  the  value  of  n '  is  much  larger  than  for  guidance  law  D 
over  most  of  the  trajectory,  as  indicated  in  Fig.  A. 6-4.  This  is  explained  by 
the  fact  that  between  seeker  measurements,  guidance  law  E  attempts  to 
completely  null  the  predicted  miss  distance  before  the  next  seeker  measure¬ 
ment  is  processed.  By  contrast,  the  other  guidance  laws  effectively  only 
attempt  to  reduce  the  predicted  miss  distance  at  each  stage.  This  law  is 
useful  in  the  sense  that  it  provides  a  lower  bound  on  the  miss  distance  that 
can  be  achieved,  for  the  assumed  airframe,  seeker  noise,  etc. 


The  computational  requirements  of  guidance  laws  B,  C,  D,  and 
E  are  primarily  dictated  by  the  Kalman  filter.  In  Section  A.  7  the  follow¬ 
ing  conservative  estimates  of  storage  capacity  and  processing  time  per 
measurement  cycle  are  obtained  for  a  typical  existing  digital  mini¬ 
computer: 
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Storage  2:  600  words 
Processing  Time  2:  0.01  sec 

The  storage  requirement  is  well  within  the  memory  capacity  of  modern 
lightweight  computers,  and  the  processing  time  is  sufficiently  small  to 
achieve  good  missile  control.  Therefore,  these  guidance  laws  are  judged 
to  be  potentially  suitable  for  mechanization  in  tactical  missiles,  provided 
their  performance  is  sufficiently  superior  to  conventional  laws  to  justify 
the  cost  of  the  computer.  The  issue  of  performance  is  treated  in  Chapter  3. 
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3.  GUIDANCE  SYSTEM  EVALUATION 


3.1  EVALUATION  PROCEDURE 

The  Covariance  Analysis  Describing  Function  Technique 
(CADET)  is  a  powerful  new  approximate  method  for  analyzing  the  statis¬ 
tical  behavior  of  nonlinear  stochastic  systems,  particularly  for  develop¬ 
ing  performance  sensitivity  curves  (Ref.  6).  Its  advantage  lies  in  the 
fact  that  mean  square  values  of  system  state  variables  for  a  given  system 
design  can  be  determined  from  one  solution  of  a  matrix  differential  (or 
difference)  equation.  By  contrast, many  solutions  of  the  system  equa¬ 
tions  of  motion  are  required  to  analyze  statistics  by  monte  carlo  tech¬ 
niques. 


The  guidance  system  model  displayed  in  Fig.  2.1-2  has  a 
mixed  continuous -discrete  character  in  che  sense  that  the  optimal 
guidance  laws  discussed  in  Chapter  2  process  data  at  discrete  times 
whereas  the  seeker,  missile,  and  airframe  equations  of  motion  are 
described  by  differential  equations.  Appendix  B  presents  a  detailed  ex¬ 
planation  of  how  CADET  is  used  to  analyze  such  a  system;  the  principal 
steps  to  be  followed  are  summarized  below: 

•  Replace  each  nonlinear  element  by  its  correspond¬ 
ing  random  input  describing  function  gain,  based 
upon  an  assumed  probability  density  function  for 
the  input  to  the  nonlinearity. 

•  Using  the  resulting  linear  system  model,  employ 
conventional  covariance  analysis  techniques  to 
propagate  the  statistics  of  the  system  state  vec¬ 
tor—  i.e.,  its  mean  and  covariance— recognizing 
that  the  describing  function  gains  are  functions  of 
those  statistics. 
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•  Compute  the  rms  miss  distance  at  the  intercept 
time  from  the  elements  of  the  system  covariance 
matrix. 


Hal 


*3 


I 

I 

l 
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Typically,  the  input  to  each  nonlinearity  is  a  state  variable  that 
is  assumed  to  be  a  gaussian  random  process.  The  gaussian  assumption 
is  based  upon  the  fact  that  the  system  state  variables  tend  to  be  a  super¬ 
position  of  past  values  of  both  the  nonlinearity  outputs  and  various  random 
system  inputs.  Thus,  in  the  sense  of  the  ’’central  limit  theorem”  (Ref. 
8),  the  probability  density  function  for  the  system  state  variables  tends 
toward  the  gaussian  form,  regardless  of  the  densities  of  individual  random 
events.  Of  course,  this  assertion  is  only  approximately  valid;  the  extent 
to  which  it  holds  depends  upon  the  amount  of  low  pass  filtering  in  the  sys¬ 
tem,  the  number  of  nonlinearities,  and  the  bandwidth  of  the  random  inputs. 
To  avoid  actually  calculating  the  required  probability  density  function,  a 
task  at  least  as  time-consuming  as  performing  monte  carlo  simulations, 
the  gaussian  assumption  stated  above  is  imposed  at  tfte  outset.  Its  validity 
for  a  specific  system  is  investigated  through  comparisons  with  selected 
monte  carlo  results. 

The  two  nonlinearities  associated  with  the  guidance  system 
investigated  in  this  study  are  the  saturation  and  sine  functions  shown  in 
Fig.  2.1-2.  The  corresponding  describing  function  gains  are  computed 
as  functions  of  their  input  statistics  in  Section  B.5. 

In  all  the  cases  treated  here,  the  mean  value  of  the  system 
state  vector  is  zero;  hence  we  are  concerned  only  with  the  propagation 
of  the  system  covariance  matrix  according  to  the  procedure  described  in 
Section  B.4,  and  calculating  the  rms  miss  distance.  The  latter  is  given 
by  the  rms  value  of  y^  (Fig.  2.1-3)  at  the  intercept  time.  The  next  sec¬ 
tion  gives  comparisons  of  guidance  law  performance  for  a  range  of 
representative  system  parameter  values. 


1 
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3.2  SIMULATION  RESULTS 

In  order  to  compare  the  performance  of  the  guidance  laws 
listed  in  Table  2. 2-1,  nominal  values  were  chosen  for  the  missile-target 
engagement  initial  conditions,  and  for  the  parameters  of  the  model  in 
Fig.  2.1-2.  These  quantities  are  listed  in  Table  3.  2-1.  All  exceptions 
to  the  nominal  conditions  are  explicitly  stated  in  the  subsequent  dis¬ 
cussion. 

Observe  that  the  autopilot  parameters  are  chosen  to  yield 
first-order  dynamics  with  a  0.1  sec  time  constant  (10  sec“*  bandwidth). 
The  seeker  noise  levels  selected  represent  an  angular  measurement 
error  of  about  0.7  mrad  at  a  missile-target  separation  of  10, 000  feet. 
The  missile-target  initial  conditions  yield  a  nearly  head-on  intercept 
trajectory. 

To  determine  the  validity  of  CADET  as  an  analysis  technique, 
the  rms  miss  distance  was  computed  using  both  monte  carlo  and  CADET 
techniques  for  selected  trajectories.  Results  are  shown  in  Fig.  3. 2-1 
for  conventional  proportional  guidance  (Law  A),  as  a  function  of  the  mis¬ 
sile  acceleration  limit;  other  off-nominal  conditions  are  stated  in  the 
figure.  The  monte  carlo  points  are  obtained  from  200  missile  trajectory 
simulations.  The  dashed  line  represents  the  value  of  rms  miss  which 
would  be  obtained  from  a  linear  covariance  analysis  that  neglects  all  the 
nonlinear  effects. 

Evidently,  the  discrepancy  between  CADET  and  monte  carlo 
results  tends  to  increase  as  amax  decreases.  This  is  to  be  expected 
because  the  system  nonlinear  behavior  becomes  more  pronounced  as 
amax  ^ecreases  CADET,  being  an  approximate  method  for  analyzing 
nonlinear  systems,  will  tend  to  have  larger  errors.  However,  the  fact 
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TABLE  3. 2-1 
NOMINAL  CONDITIONS 


Quantity 


Autopilot  Parameters  1  b. 


Target  Maneuver  Bandwidth,  xt 
rms  Target  Acceleration,  o t 
Missile  Acceleration  Limit,  a 

max 

Seeker  Time  Constant,  Tj 
Radome  Slope  Parameter,  kr 

Guidance  Navigation  Ratio  (Laws  A,  B,  and  C  only),  n ' 

Launch  Range 

rms  Heading  Error 

Closing  Velocity  (assumed  constant) 

Target  Velocity 
Interceptor  Velocity 

Receiver,  o-  /r 
r* 

Noise  Parameters^  Scintillation,  agr 

Range  Independent,  cr 


Nominal  Value 
Specification 


0 

0 

0.1  sec 
0 
0 

-1 

2 
2 


0.2  sec 
300  ft/sec 
800  ft/sec 


0. 1  sec 
0 

3 

24, 000  ft 
0. 15  rad 
4000  ft/sec 
1000  ft/sec 
3000  ft/sec 
4  x  10"8  ft"1 

4  ft 

4  x  10“4  rad 
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MISSILE  ACCELERATION  LIMIT,  omex(ftAec2) 


Figure  3.2-1  Cadet-Monte  Carlo  Comparison 
versus  Missile  Acceleration 
Capability 

that  CADET  captures  the  major  effect  of  the  nonlinearity,  and  gives 

accurate  results  when  the  system  is  moderately  nonlinear— e.g. ,  when 
o 

a  =  200  ft/sec  in  Fig.  3. 2-1 --makes  it  a  valuable  technique  for 
max 

comparing  different  guidance  system  designs. 

The  justification  for  using  CADET  lies  in  the  time  saving 
realized  relative  to  the  monte  carlo  approach.  The  monte  carlo  points 
in  Fig.  3.2-1  each  required  720  seconds  of  computer  time;  the  correspond¬ 
ing  values  from  the  CADET  results  each  require  25  seconds  of  computer 
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time.  Thus  many  more  parameter  studies  of  a  nonlinear  guidance 
system  can  be  achieved  with  CADET,  for  a  given  computational  budget. 

All  five  guidance  laws  are  compared  in  Fig.  3.2-2  over  a 
range  of  values  for  the  target  maneuver  bandwidth.  *  In  order  that  Law 
A  (proportional  guidance)  be  fairly  compared  with  the  others,  its  filter 
time  constant,  Tg  in  Fig.  2. 2-1,  is  optimized  to  yield  the  lowest  value 
of  rms  miss  distance  for  each  value  of  Xfc.  The  dispersion  between  the 
various  laws  generally  increases  with  xt.  However,  the  major  improvement 
overconventionalproportionalguidanceisachievedthroughuseof  Law  C, 
which  adds  a  control  gain  to  account  for  target  maneuvers.  Some  addi¬ 
tional  improvement  is  achieved  by  accounting  for  missile  autopilot 
dynamics  (Law  D).  Guidance  law  E  provides  a  lower  bound  on  the 
achievable  miss  distance,  neglecting  the  effect  of  the  kinematic  sine- 
function  nonlinearities  in  Fig.  2.1-2.  (The  latter  are  included  in  the 
truth  model  but  not  in  the  design  model. ) 

Observe  that  there  is  no  appreciable  difference  between  the 
performance  of  Laws  A  and  B,  although  Law  B  should  be  superior  because 
it  employs  a  Kalman  filter  to  estimate  line-of-sight  rate.  This  is  ex¬ 
plained  by  the  fact  that  most  of  the  miss  distance  associated  with  pro¬ 
portional  guidance  law  is  caused  by  target  maneuvers,  rather  than  hom¬ 
ing  sensor  measurement  noise.  This  is  illustrated  in  Fig.  3. 2-3  for  a 
case  with  no  target  maneuver  and  larger  noise  levels.  Evidently  the 
miss  distance  achieved  with  the  Kalman  filter  is  about  half  that  achieved 
with  an  optimized  first-order  filter.  However,  this  difference  is  not 

♦Guidance  Laws  B  through  D  listed  in  Table  2. 2-1  depend  upon  knowledge 
of  the  bandwidths  of  the  target  maneuver  (Xfc)  and  the  first-order  design 

model  for  the  autopilot  (Xm).  Unless  otherwise  stated  in  the  following 

discussion,  it  is  assumed  that  these  quantities  are  equal  to  the  values 
assumed  for  the  truth  model;  i.e. ,  the  optimal  guidance  laws  are 
matched  to  the  truth  model. 
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R-79S3 


Figure  3. 2-2  Guidance  Law  Performance  versus 
Target  Maneuver  Bandwidth 


Figure  3. 2-3  Performance  of  Proportional  Guidance 
Using  a  Kalman  Filter  versus  a  First- 
Order  Filter 
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significant  at  the  nominal  noise  levels,  when  target  maneuvers  are  also 
present. 


Another  aspect  of  the  comparison  in  Fig.  3. 2-2  is  the  fact 
that  Laws  A,  B  and  C  have  constant  navigation  ratios,  whereas  those 
associated  with  D  and  E  are  time-varying.  The  performance  of  the 
first  three  laws  can  be  improved  if  their  associated  values  of  n'are 
optimized,  in  the  same  fashion  as  the  filter  time  constant  associated 
with  Law  A  is  optimized  in  Fig.  3.2-3.  This  is  demonstrated  in  Fig. 
3.2-4  for  Law  B.  Generally,  n'  should  be  larger  than  the  value  three, 
which  is  obtained  in  Section  A.  6,  assuming  no  target  acceleration  and 
neglecting  autopilot  dynamics,  to  compensate  for  target  maneuvers  and 
autopilot  dynamics. 

Figure  3.2-5  displays  guidance  law  performance  as  a  function 
of  missile  maneuver  capability*.  As  approaches  the  target  rms 
acceleration  level  (300  ft/sec  ),  Laws  C,  D,  and  E  offer  marked  improve¬ 
ment  over  proportional  guidance  (Law  B).  In  addition,  the  dispersion 
between  D  and  E  increases;  this  is  attributable  to  the  fact  that  Law  E 
explicitly  accounts  for  the  acceleration  limit  whereas  D  does  not.  Hence 
the  former  offers  greater  improvement  over  the  latter  as  acceleration 
saturation  becomes  more  pronounced;  this  effect  is  also  observed  in 
Ref.  7.  . 


The  influence  of  the  missile  autopilot  time  constant  on  miss 
distance  is  demonstrated  in  Fig.  3.2-6.  As  the  latter  gets  larger,  Laws 
D  and  E,  which  explicitly  account  for  autopilot  dynamics,  offer  signifi¬ 
cantly  better  performance  than  Laws  B  and  C.  This  is  primarily 

♦Law  A  is  omitted  from  this  and  further  guidance  law  comparisons  on 
the  basis  that  it  yields  performance  quite  close  to  Law  B  when  its  filter 
time  constant  is  optimized  (see  Fig.  3.2-2)  at  the  nominal  noise  levels. 
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Figure  3.2-4  Performance  of  Proportional  Guidance 
as  a  Function  of  the  Navigation  Ratio 


Figure  3. 2-5  Guidance  Law  Performance  versus  Missile 
Maneuver  Capability 
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AUTOPILOT  TIME  CONSTANT  (»e) 

Figure  3.2-6  Guidance  Law  Performance  versus  Missile 
Autopilot  Time  Constant 

attributed  to  the  large  values  of  n '  (see  Fig.  A.  6-4)  used  in  Laws  D  and 
E  over  portions  of  the  missile  trajectory. 

Figure  3. 2-7  compares  the  performance  of  Laws  B,  D,  and  E 
as  a  function  of  target  acceleration.  Ordinarily  one  might  expect  the 
miss  distance  to  increase  more  rapidly  with  target  acceleration  than 
these  curves  indicate.  This  is  explained  by  the  fact  that  the  target 
trajectory  turns,  while  evading  the  missile,  so  that  the  component  of  its 
lateral  acceleration  normal  to  the  line-of -sight  is  reduced.  The  greater 
the  target  acceleration,  the  more  its  trajectory  turns  and  the  greater  the 
reduction  in  its  effective  maneuvering  acceleration,  normal  to  the  LOS. 

In  the  limit  when  the  target's  acceleration  becomes  directed  along  the 
line-of-sight,  it  will  have  little  influence  upon  miss  distance.  Thus,  the 
curves  in  Fig.  3. 2-7  tend  to  level  off  at  large  maneuver  levels.  This 
behavior  is  Raptured  by  the  kinematic  sine -functions  nonlinearities  shown 
in  Fig.  2.1-2;  it  would  not  be  visible  in  apurely  linear  system  analysis. 
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Figure  3.2-7  Guidance  Law  Performance  versus 
Target  Maneuver  Level 

In  all  the  cases  treated  thus  far,  the  design  of  guidance  Laws 
D  and  E  is  matched  to  the  truth  model  representation  of  missile  dynamics, 
except  for  the  kinematic  nonlinearities  mentioned  above.  An  important 
issue  remaining  to  be  investigated  is  the  sensitivity  of  guidance  law  per¬ 
formance  to  variations  in  missile  parameters  from  their  assumed  values. 

Recall  that  one  of  the  potentially  important  effects  is  autopilot- 
seeker  coupling  through  the  seeker  aberration  effect.  This  is  represent¬ 
ed  in  Fig.  2. 1-2  by  the  constant  kr,  which  has  heretofore  been  chosen  as 
zero.  Figure  3.2-8  illustrates  the  effect  on  guidance  system  performance 

25 
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Figure  3. 2-8  Guidance  Law  Sensitivity  to  Aberration  Error; 
Dominant  Autopilot  Pole  =-1.0  sec-1 


when  kr  is  nonzero.  For  this  case,  the  parameters  Oj,  a2,  bp  bg,  and 
bg  are  selected  to  yield  autopilot  poles  at  -1  rad/sec  and  -7.5  ± j  15  rad/ 
sec,  and  zeros  at  ±  20  rad/sec.  The  design  model  for  Laws  D  and  E 
assumes  the  autopilot  is  first-order  with  time  constant  (r  )  equal  to  1 
sec.  The  value  of  a^  is  chosen  to  be  representative  of  high  altitudes, 
where  autopilot  seeker  coupling  is  most  pronounced.  With  the  exception 
that  here  k  /  0,  the  above  conditions  are  nearly  the  same  as  those  in 
Fig.  3. 2-6  when  the  autopilot  time  constant  =  1.0  sec.  Figure  3. 2-9 
displays  the  sensitivity  to  kf  for  another  case  where  the  dominant  auto¬ 
pilot  pole  is  3.3  rad/sec. 

As  kr  increases,  the  performance  advantage  of  Laws  D  and  E 
deteriorates  relative  to  Laws  B  and  C.  This  is  attributed  to  the  fact  that 
the  high  value  of  n f  associated  with  Laws  D  and  E  is  incompatible  with 
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Figure  3.2-9  Guidance  Law  Sensitivity  to  Aberration  Error  ; 

Dominant  Autopilot  Pole  =  -3.3  rad/sec 

the  parasitic  attitude  loop  in  Fig.  2.1-2,  because  the  latter  is  not 
accounted  for  in  the  derivation  of  the  guidance  law.  Generally  kr  may 
be  an  unknown,  even  time -varying  quantity,  so  that  it  cannot  be  accounted 
for  exactly.  Consequently,  some  degradation  in  the  performance  of  the 
Mhigh-gainM  guidance  laws  will  be  experienced  in  situations  where  the 
effect  of  kr  is  important— i.  e. ,  at  high  altitudes.  The  performance  of 
Law  E  tends  to  degrade  more  rapidly  than  the  other  guidance  laws 
because  of  its  relatively  large  navigation  ratio  (see  Fig.  A.  6-4).  It  is 
the  most  sensitive  to  errors  in  the  design  assumptions. 

Finally,  the  sensitivity  of  guidance  law  performance  to  an  error 
in  knowledge  of  the  dominant  autopilot  time  constant  is  shown  in  Fig. 
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3. 2-10  for  Laws  D  and  E.  The  actual  time  constant  is  1  sec,  corres¬ 
ponding  to  the  conditions  at  the  right  hand  side  of  the  graph  in  Fig.  3.2-6; 
the  design  value  of  the  autopilot  time  constant  (r  )  is  varied  over  the 
range  indicated  on  the  abscissa  in  Fig.  3.2-10.  When  the  design  value  of 
rm  is  one,  both  laws  E  and  D  are  matched  to  the  autopilot  and  law  E  has 
the  best  performance.  When  rm  4-  1,  Law  E  degrades  more  rapidly  than 
Law  D  because  of  its  larger  navigation  ratio. 


DESIGN  VALUE  OF  AUTOPILOT  TIME  CONSTANT, Tm  (sec) 


Figure  3. 2-10  Guidance  Law  Sensitivity  to  Design  Value 

of  the  Autopilot  Time  Constant 
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4.  SUMMARY  AND  CONCLUSIONS 

4.1  SUMMARY 

This  report  presents  comparisons  of  performance  for  the 
guidance  laws  summarized  in  Table  4.1-1  and  derived  in  Appendix  A. 
Each  law  is  evaluated  in  a  guidance  system  model  that  incorporates  sig¬ 
nificant  linear  and  nonlinear  dynamic  effects.  These  include  the  follow¬ 
ing  missile  characteristics:  homing  seeker  noise*  homing  seeker 
dynamics,  autopilot  dynamics,  maneuvering  acceleration  limit,  parasitic 
autopilot-seeker  dynamic  coupling,  random  target  accelerations,  and 
nonlinear  missile -target  geometric  effects.  The  rms  miss  distance 
achieved  with  each  guidance  law  is  calculated  using  a  recently  developed 
analytical  procedure  for  determining  the  statistical  properties  of  non¬ 
linear  systems—the  Covariance  Analysis  Describing  Function  Technique 
(CADET).  The  latter  achieves  a  good  approximation  to  monte  carlo 
results,  and  offers  the  advantage  of  requiring  significantly  less  compu¬ 
tation  time. 

TABLE  4.1-1 

MISSILE  GUIDANCE  LAWS 


Guidance 

Laws 

Distinguishing  Characteristics 

Type  of 
Noise  Filter 

A 

Proportional  Guidance 

First-Order, 

Low-Pass 

B 

Proportional  Guidance 

Kalman 

C 

Optimal  Linear  Guidance; 

Accounts  for  Target  Maneuvers 

Kalman 

D 

Optimal  Linear  Guidance; 

Accounts  for  Target  Maneuvers 
and  Missile  Autopilot  Dynamics 

Kalman 

E 

Optimal  Nonlinear  Guidance; 
Accounts  for  Target  Maneuvers 
Missile  Autopilot  Dynamics, 
and  Missile  Airframe  Saturation 

Kalman 
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4. 2  CONCLUSIONS 

The  principal  conclusions  of  this  work  are  as  follows: 

•  Conventional  proportional  guidance  using  a  first 
order  filter  (Law  A)  is  adequate  for  nonmaneuver¬ 
ing  targets  and  moderate  homing  sensor  noise 
levels.  The  miss  caused  by  large  noise  levels 
can  be  reduced  by  as  much  as  50  percent  using 

a  Kalman  filter  (Law  B)  to  process  seeker  meas¬ 
urements  . 

•  Guidance  LawC  generally  offers  the  biggest 
relative  improvement  over  proportional  guidance 
against  maneuvering  targets,  provided  the  mis¬ 
sile  autopilot  time  constant  is  small  (<0.2  sec). 

•  Guidance  Laws  Dand  E  are  significantly  better 
than  LawC  when  the  missile  autopilot  time  con¬ 
stant  is  large  (>  0.2  sec). 

•  Laws  D  and  E  are  sensitive  to  errors  in  the 
design  assumptions  about  the  missile  autopilot 
dynamics;  LawE  is  most  sensitive  because  of  its 
relatively  large  control  gains. 

•  The  principal  computational  burden  imposed  by 
the  optimal  guidance  laws  is  associated  with  the 
Kalman  filter.  The  storage  capacity  and  meas¬ 
urement  processing  time  needed  to  mechanize 
the  filters  are  conservatively  estimated  at  600 
words  and  0.01  sec,  respectively.  These  require¬ 
ments  are  compatible  with  modem  airborne-type 
computers. 

•  CADET  is  found  to  be  a  useful  technique  in 
analyzing  guidance  system  performance.  The 
saving  in  computer  time  over  200  run  monte 
carlo  analysis  is  about  30:1  for  the  system 
model  used  in  these  studies. 

The  above  conclusions,  together  with  the  quantitative  performance 
results  provided  in  Chapter  3  will  provide  a  basis  for  choosing  among 
optimal  guidance  laws  in  particular  applications. 
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APPENDIX  A 

TRUTH  MODEL  DEVELOPMENT 


A  meaningful  comparison  of  homing  guidance  systems  for  tacr 
tical  missiles  requires  realistic  models  for  the  missile  and  its  target 
engagement  geometry  model  in  order  to  accurately  evaluate  terminal  miss 
distance.  This  model  should  include  the  important  dynamics  and, system 
nonlinearities  which  influence  performance,  and  yet  be  representative 
of  tactical  missiles  in  general.  A  system  truth  model  that  satisfies  these 
requirements  is  developed  in  this  appendix.  The  various  modules  dis¬ 
cussed  in  subsequent  sections  are  based  on  the  functional  diagram  in  Fig. 
2.1-1.  It  should  be  noted  at  the  outset  that  the  model  developed  herein 
assumes  that  the  target  and  missile  motion  are  constrained  to  a  plane. 
Consequently,  development  of  the  missile  and  guidance  models  is  limited 
to  a  single  channel. 

A.  1  SEEKER  MODULE 

The  function  of  the  seeker  subsystem  is  two -fold;  it  provides 
the  measurements  of  target  motion  required  to  mechanize  the  guidance 
law,  and  it  tracks  the  target  with  the  antenna  or  energy  receiving  device.* 
The  typical  seeker  hardware  consists  of  two  or  three  gimbals  on  which 
are  mounted  gyros  and  an  antenna.  Associated  with  each  gimbal  is  a 
servomechanism  which  is  used  to  adjust  its  angular  orientation  in  res¬ 
ponse  to  the  tracking  error  signal  measured  by  the  radar  receiver. 

(Only  one  gimbal  and  its  associated  dynamics  is  required  for  the  planar 
motion  model. )  It  should  be  noted  that  there  are  also  body  mounted 
antenna  systems  which  do  not  use  moveable  gimbals  to  position  the  anten¬ 
na.  These  systems  use  either  a  fixed  antenna  position  relative  to  the 
missile  or  electronic  beam  steering  by  means  of  a  phased  array  (radar) 
antenna.  These  configurations  are  rather  atypical  of  tactical  missiles 
and  will  not  be  specifically  considered,  although  the  use  of  electronic 

*  Generally  we  shall  use  the  term  "antenna"  to  refer  to  any  type  of  energy 
collecting  device— radar,  infrared,  or  optical. 
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beam  steering  is  in  many  ways  analogous  to  the  gimballed  system  as  far 
as  the  resulting  guidance  system  operation  is  concerned. 

Seeker  Measurement  Geometry -The  fundamental  measure¬ 
ment  obtained  from  the  homing  sensor  receiver  is  assumed  to  be  the 
indicated  angular  position  of  the  target  relative  to  the  antenna  center-line 
or  boresight.  The  guidance  laws  considered  in  this  report  require  line- 
of -sight  (LOS)  angle  or  LOS  angular  rate  as  the  fundamental  measurement 
for  terminal  guidance.  Illustrated  in  Fig.  A.  1-1,  the  LOS  angle,  0,  is 
the  angle  between  a  line  from  the  center  of  the  seeker  antenna  to  the 
target,  and  some  arbitrary  non-rotating  reference  line.  It  is  convenient 
to  select  this  reference  equal  to  the  LOS  position  at  the  beginning  of  the 
homing  guidance  phase.  Consequently,  0 (t)  at  time  t  is  the  total  change 
in  the  angular  position  of  the  LOS  relative  to  the  initial  LOS. 

The  angular  position  of  the  missile  body  center  line,  0  ,  is 
measured  relative  to  the  initial  LOS  as  shown  in  Figure  A.  1-1.  The 
angular  position  of  the  antenna  centerline  measured  relative  to  the  body 
centerline  is  defined  as  0h  in  Figure  A.  1-1.  Therefore,  0  is  given  by 

0  *  0™  +  +  e 

m  n 

Alternatively,  by  writing  Eq.  (A.  1-1)  as 

c  -  0  -  0  -  6. 

m  h 

we  obtain  an  expression  for  the  true  boresight  error.  It  is  important  to 

note  that  boresight  error  is  a  function  of  both  the  missile  attitude  relative 

to  inertial  space  and  the  angular  position  of  the  antenna  relative  to  the 

•  » 

missile  center-line.  Since  0  or  0  is  the  desired  measurement  for 
guidance  purposes,  it  is  necessary  to  remove  missile  motion  from  the 
LOS  measurement  data. 


(A. 1-1) 

(A.  1-2) 
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One  requirement  on  the  seeker  system  is  to  keep  the  antenna 
pointed  at  the  target  so  that  c  is  always  much  smaller  than  the  beam 
width  of  the  received  energy.  In  the  region  of  small  e,  the  seeker  re¬ 
ceiver  measurement  of  indicated  boresight  error  is  nearly  linear.  How¬ 
ever,  if  c  cannot  be  considered  small  relative  to  the  antenna  beam  width, 
the  receiver  boresight  error  processor  operation  may  become  nonlinear,  as 
illustrated  in  Fig.  A.  1-2.  In  fact,  if  c  is  allowed  to  approach  the  half 
beamwidth  of  the  antenna,  the  receiver  detection  eiri  uitry  will  at  some 
point  loose  lock  and  all  guidance  information  will  be  lost.  Therefore, 
the  seeker  must  track  the  target  sufficiently  closely  so  that  large  bore¬ 
sight  errors  do  not  occur.  Otherwise  the  nonlinearity  of  the  boresight 
error  position  should  be  considered  as  an  important  system  nonlinearity. 
Since  the  actual  form  of  the  boresight  error  processor  nonlinearity  is 
strongly  dependent  on  the  specific  beam  width,  processing  scheme  (mono¬ 
pulse  radar,  c.w.  radar,  etc.)  and  detector  characteristics  of  indivi¬ 
dual  systems,  it  will  not  be  included  in  the  general  system  model  for  the 
current  study.  It  will  be  assumed  that  the  beam  width  and  the  tracking 
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Figure  A.  1-2  Signal  Processing  Nonlinearities 


response  of  the  seeker  are  adequate  to  keep  the  boresight  error  pro¬ 
cessor  in  its  linear  region. 


Aberration  Error- The  aberration  angle  error  is  the  result  of 
nonlinear  distortions  in  the  received  energy  as  it  passes  through  the  pro¬ 
tective  covering  (radome  in  the  case  of  a  radar  homing  sensor)  over  the 
antenna.  This  distortion  produces  a  false  boresight  error  signal,  c ', 
which  is  interpreted  as  an  error  in  the  angular  position  and  motion  of  the 
target  by  the  guidance  package.  Referring  to  Fig.  A.  1-3,  the  indicated 


boresight  error  in  the  presence  of  aberration  angle,  @r,  is  given  by 
by 


e'  =  e  +  e  -  e  -  ek  (a.i-3) 

r  m  n 

The  size  of  this  measurement  error,  $r,  depends  on  the  orientation  of  the 
antenna  with  respect  to  the  antenna  cover,  which  is  fixed  to  the  missile 
airframe.  This  dependency  of  er  on  couples  body  motion  into  the  bore¬ 
sight  error  signal,  thus  forming  what  is  called  the  "parasitic  attitude  loop”. 
The  latter  can  drastically  alter  missile  response  characteristics  and  in 
turn  increase  miss  distance.  This  is  particularly  true  at  high  altitudes 
where  the  missile  body  motion  tends  to  be  greatest.  (This  effect  is  dis¬ 
cussed  more  fully  in  Section  A.  3). 
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Figure  A.  1-3  Effect  of  Aberration  Error 


The  aberration  angle  error  can  be  a  nonlinear  function  of 
several  factors:  the  angle  between  the  missile  center  line  and  the  LOS  to 
the  target  (look  angle)  as  illustrated  in  Fig.  A.  1-4;  the  thickness  dis¬ 
tribution,  material,  shape,  and  optical  or  electrical  properties  of  the 
antenna  cover;  frequency,  and  polarization  of  the  received  signal;  manu¬ 
facturing  tolerances;  erosion  of  the  surface  during  flight.  Therefore,  an 
accurate  model  may  require  a  nonlinear,  time-varying  statistical  charac 
terization  of  the  radome.  However,  since  these  characteristics  tend  to 
vary  over  rather  wide  limits  depending  on  the  particular  application  and 
missile  configuration,  a  constant  aberration  error  slope  model  is  used 
to  capture  the  important  body  coupling  effect,  consistent  with  the  desire 
for  a  general  system  model. 
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Figure  A.  1-4  Typical  Aberration  Angle  Error  as  a 
Function  of  Look  Angle 

A  linear  model  for  the  general  aberration  angle  characteristic 
given  in  Fig.  A.  1-4  is  obtained  from  a  simple  Taylor  series  approxima¬ 
tion 

er  *  8rb  +  ®  -  em>  kr  (A-1-4* 

where  kr  is  the  error  slope  and  erb  is  a  bias  angle.  Substituting  Eq. 

(A,  1-4)  into  Eq.  (A.  1-3)  for  c '  yields 

c '  a  (i  +  kr)  (e  -  em)  +  erb  -  eh  (A.  1-5) 

The  boresight  error  bias  will  be  assumed  negligible  relative  to  other  sys¬ 
tem  errors.  There  also  is  a  possible  contribution  of  the  aberration  error 
to  measurement  noise  when  the  frequency  of  the  received  signal  is  varied  in 
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a  pseudo -random  fashion  to  reduce  the  effect  of  a  potential  enemy  jammer, 
e.  g. ,  when  the  seeker  is  an  active  or  semiactive  radar.  This  noise  is 
simply  treated  as  a  contributor  to  range  independent  noise  (Section  A.  2). 

Seeker  Track  and  Stabilization  Loops  -  The  assumed  configura¬ 
tion  for  the  seeker  is  illustrated  in  Fig.  A.  1-5.  The  indicated  boresight 
error  ,  *  ',  is  scaled  by  1/r-p which  forms  the  desired  rate  command  to 
the  stabilization  loop.  (Although  an  actual  system  requires  the  imple¬ 
mentation  of  two  or  three  channels  to  account  for  motion  in  three  dimen¬ 
sions,  only  one  channel  is  required  for  the  planar  intercept  model  con¬ 
sidered  herein. ) 


R-9210 
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Figure  A.  1-5  Seeker  Subsystem  Block  Diagram 
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The  stabilization  dynamics  are  comprised  of  the  gimbal  servo 
and  rate  gyro  (mounted  on  the  antenna)  and  typically  have  very  wide  band¬ 
width;  e.g. ,  greater  than  100  rad/sec.  The  track  loop  model  has  simple 
first-order  dynamics;  it  commands  a  gimbal  rate  proportional  to  the 
measured  boresight  error.  The  loop  attempts  to  drive  the  boresight 
error  to  zero,  thereby  causing  the  antenna  to  track  the  target.  It  is 
straightforward  to  show  that  the  linear  transfer  function  from  8  to  c  '  — 
assuming  unity  gain  for  the  antenna  cover,  the  signal  processor,  and  the 
stabilization  dynamics  —  is 


- —  (A.  1-6) 

1  +  sr^ 

Therefore,  at  low  frequencies  (to  <  1/tj),  the  indicated  boresight  error 
is  proportional  to  the  LOS  rate.  The  latter  is  the  desired  measurement  for 
classical  proportional  navigation  guidance, which  commands  a  missile 
lateral  acceleration  proportional  to  the  LOS  rate. 

It  was  assumed  in  the  above  development  of  the  seeker  opera¬ 
tion  that  the  boresight  error  processor  nonlinearity  illustrated  in  Fig. 

A.  1-2  could  be  neglected.  Equation  (A.  1-6)  provides  an  indication  of  the 
important  region  of  boresight  error  linearity.  Using  the  fact  that  e '  is 
proportional  to  8  in  steady  state  for  a  constant  e ,  we  obtain  the  follow¬ 
ing  expression  for  c  'max: 


€ 


/ 

max 


v  i  0 


1  max 


(A. 1-7) 
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If  is  sufficiently  small,  can  be  held  within  the  linear  range  of 
the  received  beamwidth.  The  resulting  seeker  block  diagram  with  the 
linear  aberration  error  model  is  given  in  Fig.  A.  1-6. 

The  importance  of  the  aberration  angle  error  slope,  kf,  is 
illustrated  by  the  linear,  continuous  transfer  function  relating  0^  to  0 ' 
in  Fig.  A.  1-6.  The  latter  is  given  by 


-k 


(1  +r^s) 


(A.  1-8) 


Thus  the  measured  LOS  rate  is  corrupted  by  a  term  proportional  to  body 
rate.  Since  body  rate  is  a  result  of  commanded  acceleration,  a  loop  is 
formed  which  can  have  a  destabilizing  effect  on  missile  attitude  and  in¬ 
crease  miss  distance.  Note  that  when  kr  is  zero,  the  contributions  from 
the  body  angular  rate  input  in  Fig.  A.  1-6  cancer,  producing  no  effect  on  € 


R-7152o 
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ERROR  BORESIGHT  TIME  CONSTANT 
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0* 
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Figure  A.  1-6  Seeker  Model  with  Track  Loop 
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LOS  Angle  Measurement  —  Although  classical  proportional 
guidance  uses  measurements  of  LOS  rate,  it  is  more  convenient  to  use 
measurements  of  LOS  angle  in  guidance  laws  that  utilize  a  Kalman  filter. 
To  obtain  the  latter,  we  define  the  measured  LOS  angle,  0 ',  as 

0'  =  (1  +kr)  9-kr0m  (A.  1-9) 


Then  from  Fig,  A.  1-6  it  follows  that 

T.  S  0  / 


(A. 1-10) 


Since  the  boresight  error  is  an  observable  quantity,  Eq.  (A.  1-10)  can  be 
inverted  to  yield 


1  +ST. 

e'«— (a.  l-ii) 

Tl® 

That  is,  0'  can  be  recovered  from  the  measured  seeker  boresight  error. 
This  is  represented  in  Fig.  2,1-2  by  direct  pickoff  of  6'. 

A.  2  NOISE  MODULE 

As  indicated  in  Section  A.  1,  LOS  angle  is  the  fundamental 
quantity  measured  by  the  Seeker  Module.  These  measurements  will 
generally  be  cornqrted  by  various  types  of  noise  which  can  be  categorized 
according  to  the  dependency  of  their  rms  levels  on  the  missile-to-target 
range.  The  actual  noise  levels  and  bandwidths  are  dependent  on  the 
exact  form  of  the  measurement  signal  processor,  target  configuration 
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and  characteristics,  environmental  conditions  and  a  multitude  of  related 
system  effects.  However,  using  measurements  obtained  from  actual 
hardware  or  mathematical  models,  most  of  the  observed  measurement 
noise  can  be  lumped  into  one  of  three  assumed  forms:  receiver*  noise, 
range  independent  noise  or  angular  scintillation  noise. 

Receiver  Noise  -  Receiver  noise  consists  primarily  of  thermal 
noise  generated  by  the  antenna  and  receiver  electronics  on  board  the 
missile.  The  effective  amplitude  of  this  noise  increases  with  increasing 
range  because  of  the  corresponding  decreasing  signal-to -noise  ratio. 
There  are  in  general  three  types  of  missile  receivers  which  can  be  con¬ 
sidered: 


•  Passive  —  Target  supplies  radiated  signal 

•  Semi-Active  —  Target  is  illuminated  by  a  source 

which  is  not  on  board  the  missile 

•  Active  —  Target  is  illuminated  by  source  on 

board  the  missile 

The  receiver  will  generally  include  some  type  of  automatic 
gain  control  which  attempts  to  keep  the  receiver  signal  power  nearly 
constant.  As  a  result,  the  effective  noise  level  will  change  with  received 
signal  power  relative  to  some  reference  level.  A  normalized  angular 
measurement  noise  model  will  be  defined  which  uses  the  variance  (or 
power  spectral  density)  of  the  indicated  boresight  error,  measured  at  a 
range  that  yields  a  signal-to-noise  ratio  of  unity  as  the  reference  level. 
The  resulting  expressions  for  the  variances  of  the  effective  noise  on  the 
LOS  angle  for  the  three  types  of  receivers  are: 

♦Again,  the  term  receiver  can  refer  to  any  type  of  homing  sensor. 
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Passive 


2  2  /  rtm 

a  =  a  I  - 

r  u  rp  ^  rQ 


(A.  2-1) 


Semi-Active 


a  =a 
r  rsa 


(T  tmrit\ 


(A.  2-2) 


Active 


where: 


or  =  a 

r  ra 


(A.  2-3) 


rtm  =  missile"to“target  range 
r.^  =  illuminator-to  -target  range 

Tq  =  range  at  which  S/N  =  1 
2 

CTrp  =  P*83*76  noise  variance  at  S/N  =  1 
2 

ct  „  =  semi-active  noise  variance  at  S/N  =  1 
rsa 

2 

a  „  =  active  noise  variance  at  S/N  =  1 
ra 

cr^.  =  effective  target  crossection 
Oq  =  reference  target  crossection 

All  three  types  of  receiver  noise  exhibit  the  characteristic  of  increasing 
variance  with  increasing  range.  Also,  note  that  the  passive  model  is  not 
dependent  on  the  effective  target  crossection  since  the  energy  received 
by  the  missile  is  direct  radiation  and  not  the  result  of  reflected  energy. 
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Often  the  illuminator  to  target  range  for  a  semi-aetive  radar 
is  nearly  constant;  in  this  case  the  semi-active  noise  variance  will  have 
nearly  the  same  variation  with  r^  as  the  passive  noise  for  a  given 
target  crossection.  Since  most  tactical  missiles  tend  to  be  passive  or 
semi-active,  the  passive  noise  model  will  be  used  in  the  truth  model. 

The  noise  bandwidth  is  dictated  by  the  post-detection  bandwidth 
of  the  receiver  which  in  general  tends  to  be  much  larger  than  the  signal 
bandwidth.  Consequently,  it  can  be  assumed  that  the  noise  is  "white" 
over  the  signal  spectrum  of  interest  without  loss  of  generality.  In  the 
case  of  a  sampled  data  system,  errors  in  the  sampled  and  held  values  of 
LOS  angle  are  assumed  to  be  uncorrelated  from  sample  to  sample.  If  it 
is  assumed  that  the  receiver  has  a  double-sided  square  spectrum  Af  Hz 
in  width  with  a  constant  noise  power  spectral  density  of  qr  rad^/Hz, 
the  variance  of  the  noise  is  simply 

a2  =  qr  Af  (rad2)  (A.  2-4) 

If  the  receiver  output  is  sampled  and  held  with  a  sample  period 
of  rg  sec,  where  Tg  is  much  larger  than  1/Af,  resulting  in  independent 
noise  samples,  the  output  noise  spectral  density,  qg  ,  has  a  low  frequency 
level  of  approximately 

<is  *  Or  Af  Ts  =  °r  Ts  (tH — )  (A-  2-5> 

Assuming  that  the  bandwidth  of  the  signal  is  smaller  than  the  sampling 
frequency,  Eq.  (A.  2-5)  yields  the  level  of  the  equivalent  white  noise 
spectral  density. 

Range  Independent  Noise  -  Range  independent  noise  is  a  collec¬ 
tion  of  all  noise  sources  which  contribute  a  constant  rms  error  in  the 
measurement  of  LOS  angle  throughout  the  flight.  Typical  sources  include 
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servo  noise  generated  by  the  seeker  servo,  and  target  amplitude  scintilla¬ 
tion  noise.  It  will  be  assumed  that  the  noise  is  also  white  over  the 
receiver  bandwidth,  and  the  preceeding  discussion  of  bandwidths  and 

sampling  is  applicable.  The  variance  of  the  range  independent  noise  is 
2 

defined  as  or  . 


Angular  Scintillation  Noise  -  Angular  scintillation  noise  is 
caused  by  the  wandering  of  the  apparent  centroid  of  radiation  across  the 
visible  surface  of  the  target.  Although  the  magnitude  of  the  wander  is 
essentially  range  independent,  the  equivalent  noise  on  the  LOS  angle 
measurements  increases  as  range  decreases,  resulting  in  a  measurement 
noise  variance  of 


ffs  = 


(A.  2-6) 


where 


cr^d  =  variance  of  the  apparent  wander  distance 

Total  Measurement  Noise  The  total  measurement  noise 
variance  is  the  sum  of  the  variances  of  the  individual  uncorrelafced  noise 
components, 


2  2  2,2 

ae  =  ar  +  0>  +  as 


(A.  2-7) 


As  previously  noted,  receiver  and  range  independent  noise  are  generally 
assumed  wide-band  relative  to  the  guidance  system  bandwidth.  Angular 
scintillation  noise  is  in  general  a  narrow  band  source  and  is  often  modeled 
as  white  noise  through  a  low  pass  filter  with  a  time  constant  which  depends 
primarily  on  the  target  motion  spectrum.  On  the  other  hand,  if  for  a 
radar  seeker  the  radar  frequency  is  changed  in  a  pseudo-random  manner 
from  sample-to-sample  so  as  to  reduce  the  ability  of  the  enemy  to  jam 
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the  missile  receiver,  then  the  apparent  centroid  of  target  radiation  will 
tend  to  be  independent  from  sample-to-sample.  It  will  be  assumed  that 
all  three  noise  sources  are  independent  and  wide-band  relative  to  the 
guidance  system  bandwidth. 

A.  3  AUTOPILOT  MODULE 

The  function  of  the  autopilot  subsystem  is  three-fold;  it  provides 
the  required  missile  lateral  acceleration  response  characteristics,  it  must 
stabilize  or  damp  the  bare  airframe,  and  it  reduces  the  missile  perform¬ 
ance  sensitivity  to  disturbance  inputs  over  the  required  flight  envelope. 

The  autopilot  configuration  illustrated  in  Fig.  A.  3-1  uses 
accelerometer  feedback  so  as  to  control  the  lateral  acceleration  of 
the  missile.  Lateral  acceleration  control  is  used  in  accordance  with  the 
proportional  navigation  guidance  law  which  requires  a  missile  lateral 
acceleration  proportional  to  the  measured  missile-to-target  line-of-sight 
rotation  rate.  The  body  mounted  rate  gyro  senses  body  attitude  rate, 

0  ,  which  is  used  by  the  autopilot  to  increase  the  effective  damping  ratio 
of  the  airframe  short  period  poles. 

The  aerodynamic  characteristics  of  the  missile  airframe  are 
an  integral  part  of  the  autopilot  design  and  operation.  Therefore,  the 
design  of  an  autopilot  must  be  tailored  to  each  individual  missile  airframe 
configuration  and  its  associated  aerodynamic  characteristics,  which  are 
nonlinear  functions  of  missile  velocity,  angle-of -attack,  control  surface 
deflection  and  altitude. 

It  is  standard  practice  in  the  design  of  missile  autopilots  to 
utilize  the  linearized  second  order  airframe  model  given  in  Chapter  8  of 
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Figure  A.  3-1  General  Missile  Autopilot  Configuration 


Hef.  5.  The  required  stability  derivatives  are  obtained  from  the  non¬ 
linear  moment  and  force  coefficients  by  making  the  following  assumptions: 

•  Constant  missile  velocity 

•  Body  lift  force  is  a  linear  function  of  the  change  in  angle 
of  attack,  a,  about  some  trim  condition,  a0 

•  Constant  altitude 

•  Constant  center  of  pressure 

•  Fixed  missile  mass  and  inertia 

•  Control  surface  lift  force  is  a  linear  function  of  control 
surface  deflection  angle,  6,  and  independent  of  ou 

Although  the  above  assumptions  appear  to  be  rather  restrictive,  they 
simplify  the  autopilot  design  task  considerably.  Practical  experience  has 
shown  that  the  resulting  autopilot  response  characteristics  with  the 
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nonlinear  airframe,  are  closely  approximated  by  the  linearized  response 
characteristics  near  the  given  nominal  conditions,  for  a  properly  designed 
autopilot.  In  this  report,  a  realistic  autopilot  model  is  developed  which 
requires  knowledge  of  very  few  specific  aerodynamic  parameters,  yet  its 
response  characteristics  are  easily  related  to  the  important  missile  aero¬ 
dynamic  properties. 

Acceleration  Command  Limiting  -  The  most  important  nonlinear 
characteristic  associated  with  the  airframe  is  acceleration  saturation, 
which  occurs  when  the  missile  attempts  to  pull  a  large  angle  of  attack.  It 
is  desirable  to  avoid  a  large  angle  of  attack  since  the  associated  drag 
results  in  a  rapid  loss  of  missile  velocity.  There  is  also  the  airframe 
structural  limit  which  must  not  be  exceeded.  It  is  common  practice  to 
limit  the  commanded  lateral  acceleration  so  as  to  prevent  both  angle -of- 
attack  saturation  and  structural  failure.  Therefore,  autopilot  command 
limiting  is  assumed  to  be  the  dominant  nonlinear  effect  and  all  other  non¬ 
linear  characteristics  ouch  as  actuator  angle  and  angle  rate  limiting, 
aerodynamic  nonlinearities,  instrumentation  nonlinearities,  etc.,  are 
assumed  to  be  secondary  or  equivalently  represented  as  acceleration 
limiting,  or  as  changes  in  autopilot  dynamics.  The  resulting  model  is 
simple  and  generally  applicable  to  a  wide  range  of  missile  systems,  and 
captures  what  is  known  to  be  a  dominant  nonlinear  system  characteristic 
and  an  important  factor  in  miss  distance  performance— lateral  accelera¬ 
tion  saturation. 

Linear  Autopilot  Response  Characteristics -Using  a  linearized 
airframe  model  (Ref,  5),  the  closed  loop  transfer  function  for  the  general 
autopilot  configuration  of  Fig.  A. 3-1  can  be  developed  for  specific  gains 
and  compensation,  an  example  of  which  is  given  in  Fig,  A.  3-2.  The  linear¬ 
ized  airframe  transfer  functions,  as  given  by  Eq.  (8. 1-1)  of  Ref.  5,  are: 
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Figure  A.  3-2  Linear  Autopilot  Block  Diagram 
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where 

“a  =  bare  air£rame  “tal 

£a  =  bare  airframe  damping  ratio 

a  =  developed  lateral  acceleration 

a„  =  measured  lateral  acceleration 

a 

6L  =  body  attitude  rate 

m 

6  =  control  surface  deflection 

a.  =  transfer  function  parameters  obtained  from 
linearized  aero-data  and  given  in  Eq.  (8.1-8) 
of  Ref.  5. 

The  Gg  transfer  function  in  Eq.  (A.  3-2)  will  be  assumed  equal 
to  Gp  this  is  true  if  the  accelerometer  is  located  at  the  c.  g.  of  the 
missile  and  has  a  wide  bandwidth.  The  actuator  dynamics  are  assumed 
to  be  at  a  frequency  which  is  much  higher  than  the  crossover  frequency 
of  the  rate  loop  and  will  therefore  be  ignored. 


Using  Eqs.  (A.3r2)  through  (A.  3-4)  the  closed  loop  transfer 
function  from  commanded  to  developed  lateral  acceleration  Fig.  A.  3-2 


is 


1  +  a<  s  +  a«  s2 

1  'iff 


l+bjS  +  bgS  +bgS 

'  <  2 
1  +  aj  s  +  ao  s 


(A.  3-5) 
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b3  =  — ^”5”  (A.  3-8) 

H  ®2 

k  =  closed  loop  d-e  gain  in  Fig.  A.  3-2. 


In  the  case  of  a  tail-controlled  missile,  the  transfer  function  in  Eq.  (A.  3-5) 
is  non-minimum  phase. .  .i.e. ,  it  has  a  right-half-plane  zero;  this  produces 
the  "wrong-way"  or  "tail-wags -dog"  effect,  discussed  in  Ref.  5.  The 
transfer  function  from  commanded  acceleration  to  developed  body  attitude 
rate  is: 


(1  +  a3  s) 


m  1  +  bg  s  +  bg  s  +  bg  s'* 


(A.  3-9) 


where  vm  is  the  missile  velocity. 


Recall  from  the  discussion  in  Section  A.  1  that  radome  errors 
can  couple  body  angle  rate,  6m,  into  the  LOS  measurements.  This 
forms  a  "parasitic  attitude  loop"  since  LOS  data  is  used  to  form  the 
lateral  acceleration  command  which  results  in  a  change  of  body  attitude 
in  Eq.  (A.  3-9).  As  this  coupling  becomes  large,  guidance  stability  is 
compromised  and  miss  distance  will  fend  to  increase.  An  important 
measure  of  the  required  missile  attitude  rate  is  &g ,  which  tends  to  in¬ 
crease  with  altitude  thus  intensifying  the  parasitic  attitude  loop  problem. 
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There  is  a  wide  range  of  values  that  can  be  selected  for  ojg  and 
£2  in  Eq.  (A.  3-5).  For  many  systems  it  is  possible,  and  desirable  where 
attitude  loop  stability  is  a  consideration,  to  design  the  autopilot  such  that 
a>2  is  much  larger  than  and  0*5*  The  value  of  a>2  depends  pri¬ 
marily  on  the  rate  loop  gain  for  the  assumed  configuration  of  Figure 
A.  3-2.  If  the  actuator  bandwith  is  much  larger  than  the  open  rate  loop 
crossover  frequency,  a>r,  which  is  in  turn  much  larger  than  the  airframe 
natural  frequency,  to  ,  then  to  can  be  approximated  as 

a  r 


“r  ”  1^8  m6 


where  ks  is  defined  in  Fig.  A.  3-2  and  m^  is  the  pitch-moment-effectiveness 

of  the  control  surface.  The  magnitude  of  m^  tends  to  increase  with  Mach 

number  and  decrease  with  altitude.  If  k  is  adjusted  (adaptively)  so  as  to 

s 

keep  tor  nearly  constant,  then  under  the  stated  condition,  u>2  and  will 
tend  to  remain  nearly  constant  over  a  wide  range  of  altitude  and  velocity 
conditions.  At  altitudes  less  than  10, 000  feet,  tends  to  be  nearly  con¬ 
stant  over  a  wide  range  of  missile  Mach  numbers.  At  higher  altitudes,  ajj 
tends  to  decrease  with  increasing  altitude  and  increase  with  increasing 
Mach  number.  However  for  an  assumed  constant  velocity  intercept  with¬ 
in  a  narrow  altitude  band,  will  remain  nearly  constant.  Therefore, 
the  guidance  system  truth  model  will  use  a  constant  coefficient  autopilot 
having  the  assumed  form  of  Eqs.  (A.  3-5)  and  (A.  3-9),  as  shown  in  Fig.  A.3-3. 


A.  4  GEOMETRY  MODULE 

The  general  form  of  the  intercept  geometry  is  illustrated  in 
Fig.  A.  4-1 .  The  reference  coordinate  system  is  defined  with  its  x-axis 
along  the  original  missile -to -target  line-of-sight  (LOS)  defined  at  the 
initiation  of  the  terminal  phase.  It  is  assumed  that  missile -target  motion 
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Figure  A.  3-3  Simulation  Autopilot  Model 
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Figure  A.  4-1  Intercept  Geometry 
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is  constrained  to  the  intercept  plane,  defined  by  the  orthogonal  reference 
x  and  y  axes.  The  intercept  plane  can  in  general  assume  any  orientation 
with  respect  to  the  earth,  depending  on  the  initial  relative  missile  and 
target  position  and  the  orientation  of  their  respective  velocity  vectors. 

A  nominal  co-altitude  intercept  trajectory  is  assumed  herein.  This 
eliminates  the  specification  of  an  altitude  profile  and  the  associated  re¬ 
quirement  to  handle  nonstationary  altitude  dependent  aerodynamic  charac¬ 
teristics  in  the  autopilot  model.  It  is  also  assumed  that  gravity  compen¬ 
sation  is  used  in  the  missile  guidance  law  to  negate  the  effect  of  gravity 
on  missile  performance.  Of  course,  for  any  orientation  of  the  intercept 
plane  that  is  not  horizontal,  there  will  be  an  altitude  change  for  motfon  in 
the  y-direction. 

Equations,  of  Motion  -  As  illustrated  in  Fig.  A.  4-1,  the 
reference  coordinate  system  is  non-rotating  but  is  allowed  to  translate 
with  the  missile.  This  selection  of  reference  system  is  dictated  primarily 
by  the  proportional  navigation  guidance  law  which  attempts  to  null  the 
LOS  rotation  rate,  0,  by  commanding  a  proportional  missile  acceleration 
normal  to  the  LOS.  If  this  guidance  scheme  works  well,  8  will  remain 
near  its  initial  zero  value  and  the  severity  of  a  number  of  geometric  non- 
linearities  can  be  reduced  and  in  many  cases  linearized. 

The  missile  lead  angle,  0^ ,  and  the  target  aspect  angle,  0a, 
in  Fig.  A.  4-1  define  the  orientation  of  the  respective  missile  and  target 
velocity  vectors  in  the  intercept  plane  relative  to  the  original  LOS. 
Closing  velocity,  v  ,is  defined  as  the  relative  velocity  measured  along 

w 

the  LOS;  viz., 

\  -  vm  cos  <««  +  «he  -  +  vt cos  (9a  +  ®> 


(A.  4-1) 
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where 

v  -  missile  velocity 
m 

vt  =  target  velocity 
6he  =  iiea<^nS  error 


Heading  error,  @he,  is  the  angular  error  in  the  collision-course  triangle 
defined  at  the  initiation  of  the  terminal  phase,  The  sides  of  the  collision- 
course  triangle  are  established  by  the  closing  velocity  along  the  original 
LOS,  the  missile  velocity  vector  and  the  target  velocity  vector.  Given  a 
target  aspect  angle,  0  ,  the  collision-course  missile  lead  angle  must  be 

2i 


(A. 4-2) 


If  the  orientation  and  magnitude  of  the  velocity  vectors  were  to  remain 
fixed  for  the  remainder  of  the  terminal  phase,  the  two  vehicles  would 
collide.  In  practice,  it  is  not  possible  to  achieve  the  collision  course 
lead  angle  and  the  difference  between  0^  and  the  actual  lead  angle  9^  is 
defined  as  the  heading  error.  For  missile  systems  having  a  midcourse 
phase  preceding  the  terminal  phase,  heading  error  tends  to  be  small, 
having  an  rms  value  of  a  few  degrees  or  less;  in  a  "dog-fight  missile" 
engagement,  heading  error  can  be  tens  of  degrees.  It  is  assumed  in  the 
following  model  development  that  heading  error  is  less  than  10  deg. 


The  missile  and  target  accelerations  are  assumed  to  be  normal 
to  their  respective  velocity  vectors.  In  reality,  there  are  also  components 
of  acceleration  developed  along  their  velocity  vectors  due  to  aerodynamic 
and  induced  drag  which  integrate  into  a  net  reduction  of  the  velocity  vector 
magnitudes.  This  is  particularly  true  of  the  missile  after  engine  burn-out. 
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However,  if  it  is  assumed  that  drag  acceleration  is  small  or  a  sustainer 
engine  is  used  to  maintain  a  nearly  constant  velocity,  only  the  acceleration 
normal  to  the  velocity  vector,  a  ,  need  be  considered.  Since  it  will  be 
assumed  that  the  magnitude  of  the  velocity  vector  is  constant,  the  rate  of 
change  of  its  orientation  is  proportional  to  the  normal  acceleration,  resul+  - 
ing  in  the  following  expression  for  lead  angle  rate: 


V  t 


(A.  4-3) 


The  analogous  expression  for  target  output  angle  rate  due  to  target 
acceleration,  a^.,  normal  to  its  velocity  vector  is 

6  =  ~  (A.  4-4) 

a  vt 

Therefore,  the  rotation  rates  of  the  velocity  vectors  are  proportional  to 
the  respective  lateral  accelerations. 


The  relative  velocity  normal  to  the  original  LOS,  y^,  is 


y ,  =  v,  sin  0  -  v  sin  0. 
d  t  am  l 


(A.  4-5) 


where 


)  =  /  a.dt+0  +0. 

I  v  /  &  ic  he 

m  JQ 


(A.  4-6) 


■H 


a,  dt  +  0 
t  ac 


(A.  4-7) 
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From  Fig.  A.  4-1,  8  is  given  by 


6  =  tan"1  (  (A.  4-8) 

\  xtm  / 

where  x^mis  themissile-to-target  range  measured  along  the  original  LOS. 
As  previously  mentioned,  B  tends  to  remain  small  if  the  guidance  law  is 
operating  properly  throughout  the  duration  of  the  terminal  phase.  How¬ 
ever,  during  the  last  fraction  of  a  second  as  intercept  is  approached, 
approaches  zero  and  6  will  tend  to  increase  for  a  given  value  of  y^.  With 
only  a  fraction  of  a  second  left  before  the  end  of  the  flight,  there  isn’t 
enough  time  for  the  missile  lateral  position  to  respond  to  the  rapidly  in¬ 
creasing  LOS  angle.  In  fact,  the  high  LOS  rate  generally  results  in 
saturation  of  the  autopilot  acceleration  command  before  B  becomes  large. 
Therefore,  for  modeling  purposes,  accurate  computation  of  6  is  not 
required  during  the  period  in  which  if  becomes  large,  thus  allowing  the 
small  angle  approximation  of  Eq.  (A.  4-9)  to  eliminate  the  tan  nonlinear 
operation. 

yd 

8  2=  r~-  (rad)  (A.  4-9) 

xtm 

It  is  important  to  note  that  xtm  appears  in  the  denominator  of 
Eq.  (A.  4-9).  Therefore,  if  Xjm  is  modelled  as  a  system  state,  it  is 
necessary  to  handle  the  associated  nonlinear  ratio  of  states.  However, 
if  Xj.m  can  be  modelled  as  a  deterministic  quantity  which  is  only  a  function 
of  time,  Eq.  (A.  4-9)  is  linear  and  the  computation  of  8  is  considerably 
simplified.  It  will  be  shown  in  the  following  paragraph  that  the  difficulties 
associated  with  the  computation  of  8  and  miss  distance  are  all  related  to 
a  non -constant  closing  velocity. 


A-26 


56 


THE  ANALYTIC  SCIENCES  CORPORATION 


Miss  Distance  -  Miss  distance  is  defined  as  the  point  of  closest 
approach  between  the  missile  and  target,  illustrated  in  Fig.  A.  4-2,  and 
is  given  by  tue  relative  range  when  the  closing  velocity  is  equal  to  zero 
in  Eq.  (A.  4-1).  Since  closing  velocity  is  defined  as  the  relative  velocity 
between  the  missile  and  target  projected  onto  the  instantaneous  LOS,  miss 
distance  is  equal  to  the  missile  to  target  range  when  the  inner  product  of 
range  and  range  rate  is  equal  to  zero,  viz. , 


miss 


where 


itm- '  xtm—  +  ydi 

2  =  unit  vector  in  the  x  direction 

j  =  unit  vector  in  the  y  direction 


(A. 4-10) 
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At  the  point  in  space  where  the  orthogonality  of  the  relative  range  and 
velocity  vectors  required  by  Eq.  (A.  4-10)  exists,  miss  distance  is  given 

by 


miss  =  yd 


(A. 4-11) 


For  many  intercept  situations,  x^  is  much  larger  than  y^  in  Eq.  (A.  4-11) . 
Therefore  the  downrange  component  of  miss  distance  is  neglected  in  these 
studies. 


Although  miss  distance  is  closely  approximated  by  y^,  it  (the 
miss)  is  defined  at  oniy  one  point  during  the  flight— namely  the  point 
of  closest  approach.  Therefore  it  is  conditioned  on  a  spatial  rather 
than  a  time  relationship.  In  a  monte  carlo  study,  the  orthogonality  con¬ 
dition  is  tested  near  the  expected  end  of  each  flight  and  yd  is  saved  when 
the  condition  is  satisfied.  Therefore,  if  the  closing  velocity  is  not  con¬ 
stant  or  there  is  a  significant  down  range  miss  component,  the  time-of- 
flight  at  which  miss  distance  is  measured  will  not  be  constant  from  one 
trial  to  the  next.  The  ensemble  miss  distance  statistics  are  computed 
using  the  "saved”  values  of  miss  distance,  each  of  which  occurs  at  a 
different  time  of  flight.  By  contrast,  the  covariance  analysis  technique 
provides  ensemble  performance  statistics  as  a  function  of  time.  There¬ 
fore,  additional  computations  or  approximations  are  required  to  obtain 
miss  distance  using  the  covariance  analysis  approach. 


The  primary  contributor  to  a  nonconstant  time -of -flight,  for  a 
given  initial  range,  is  the  fact  that  the  closing  velocity  in  Eq.  (A.  4-1)  is 
not  constant.  Since  it  has  already  been  assumed  that  the  magnitudes  of 
the  missile  and  target  velocity  vectors  ar  :onstant,  any  change  in  y^  is 
a  result  of  the  change  in  the  orientation  of  the  velocity  vectors.  However, 
if  the  change  in  these  angles  is  assumed  to  be  small  during  the  flight,  the 
closing  velocity  will  remain  nearly  constant  except  as  the  point  of  mini¬ 
mum  range  is  approached  where  it  rapidly  approaches  zero.  The  projec- 
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tion  of  misBlle  and  target  velocity  onto  the  original  LOS,  xfcm,  will  also 
be  nearly  constant  throughout  the  flight.  Therefore,  x^m  is  assumed  to 
be  a  deterministic  quantity  which  goes  to  zero  at  a  known  fixed  time. 

This  assumption  eliminates  the  previously  mentioned  computational  prob¬ 
lem  associated  with  6  in  Eq.  (A.  4-9). 


From  Eq.  (A.  4-10)  it  follows  that  miss  is  defined  when 


Vd 

xtm 


(A.  4-12) 


The  incremental  time,  At,  between  when  Eq.  (A.  4-12)  is  satisfied  and 
xtm  is  e^ua^  to  zero  is  approximately  Simply  dividing  Eq.  (A.  4rl2)  by 

xtm* 


A 


Yw  y 


d'd 


(A.  4-13) 


For  most  intercepts,  A  t  is  on  the  order  of  milliseconds.  The  change  in 
yd  over  such  a  small  interval  is  negligible  and  miss  distance  can  be 
approximated  as  y^  measured  when  xtm  equals  zero. 


A.  5  TARGET  MODULE 

The  target  velocity  vector  is  assumed  to  have  a  constant  mag¬ 
nitude^  and  a  direction  described  by  6  in  Fig.  (A.  4-1).  It  is  further 
assumed  that  the  target  may  have  random  changes  in  its  acceleration 
normal  to  its  velocity  vector.  The  assumed  acceleration  time  history 
model  is  a  randomly  reversing  poisson  square  wave  as  illustrated  in 

o 

Fig.  A.  5-1.  This  square  wave  switches  between  ±0  ft/sec  with  random 
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poison  distributed  switching  times  having  an  average  of  *  aero-crossings 

per  second.  The  autocorrelation  function  (Ref.  9)  for  observation  times 
taand^ls 


*(ta-tb)^V2i'|ta-tb| 

Equation  (A.  5.®  indicates  that  the  mean-squared  value  of 


(A.  5-1) 


V  approaches  zero,  at  approaches  a  constant.  The 


associated  with  afc  is 


*(oo)  =  2 


ue  of  afc  is  /r  and,  as 
power  spectral  density 


(A.  5-2) 


Figure  A.  5-1  Poisson  Square  Wave  for 


Target  Acceleration 
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where  X^  is  defined  as  the  maneuver  bandwidth  and  is  equal  to  2i/.  The 
resulting  maneuver  model  has  the  same  mean  and  autocorrelation  function 
as  a  first-order  markov  process.  Therefore,  the  target  lateral  maneuver 
acceleration  will  be  modelled  as  the  output  of  a  low  pass  filter  with  band¬ 
width  Xt  driven  by  white  noise  with  a  power  spectral  density  of  2X|?2, 
resulting  in  a  rms  acceleration  level  of  a. 

It  is  interesting  to  note  that  although  the  autocorrelation  func¬ 
tion  and  the  corresponding  power  spectral  density  for  the  poisson  square 
wave  are  identical  to  that  of  a  markov  process,  the  associated  probability 
density  functions  can  be  quite  different.  It  is  obvious  from  Fig.  A.  5-1 
that  the  poisson  square  wave  can  only  take  on  values  of  ±0.  The  resulting 
bi-polar  amplitude  probability  density  function  consists  of  impulses  with 
a  weighting  of  0. 5  at  plus  and  minus  /$  whereas  the  markov  process  is 
generally  assumed  to  have  a  gausaian  amplitude  distribution.  Therefore, 

the  response  of  an  amplitude  dependent  nonlinear  operator  could  be  quite 
different  when  driven  by  each  of  these  two  signal  forms.  However,  if  the 
random  square  wave  is  passed  through  a  narrow  band  filter  or  integrator, 
it  would  experience  rounding  due  to  the  finite  bandwidth.  In  the  case  of 
an  integrator,  the  resulting  wave  shape  would  be  a  series  of  constant 
slope  segments.  By  application  of  the  central  limit  theorm,  as  illus¬ 
trated  in  Ref.  8*,  the  resulting  output  distribution  approaches  the  gaussian 
density  function. 
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A.  6  GUIDANCE  MODULE 

The  operation  of  the  guidance  module  may  be  separated  into 
two  cascaded  functions:  1)  filtering  of  the  noisy  measurements  obtained 
from  the  seeker  and  2)  utilization  of  the  filtered  measurements  to  con¬ 
trol  the  missile  lateral  acceleration,  There  are  a  number  of  filtering 
and  control  schemes  that  can  be  used  by  tactical  missiles,  as  reported 
in  Refs.  5  and  7.  However,  the  current  study  will  be  limited  to  the  con¬ 
figurations  listed  in  Table  2. 2-1;  each  of  these  is  discussed  below. 

Configuration  A  -  This  configuration  can  be  considered  as  the 
classical  approach  to  missile  guidance.  The  boresight  error  signal  from 
the  seeker  module  is  processed  through  a  constant  bandwidth  low-pass 
noise  filter.  For  modelling  purposes,  it  is  assumed  that  the  measurement 
rate  is  much  faster  than  the  noise  filter  bandwidth.  The  resulting  system 
can  therefore  be  characterized  by  a  continuous  model  which  is  represen¬ 
tative  of  a  missile  with  either  an  analog  on-board  signal  processor  or  a 
digital  unit  with  a  high  cycle  rate.  This  system  can  also  be  considered 
as  one  of  minimum  complexity. 


The  boresight  error  signal  obtained  from  the  seeker  module  is 
a  noisy  measurement  of  LOS  rate  at  frequencies  below  that  of  the  seeker 
track-loop  bandwidth.  Therefore,  the  output  of  the  guidance  noise  filter 
is  also  a  band-limited  indication  of  LOS  rate.  Classical  proportional 
guidance  requires  the  development  of  a  lateral  acceleration  of  the  missile 
normal  to  the  LOS,  which  is  proportional  to  the  LOS  rate  as  given  by 
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where 

vc 

• 

0 

0 

n' 


closing  velocity  (ft/sec) 
LOS  rate  (rad/sec) 

LOS  angle  (rad) 
navigation  ratio 
lead  angle  (rad) 


The  proportionality  factor  is  comprised  of  the  navigation  ratio, 
closing  velocity  multiplier  and  a  geometric  gain  factor  which  accounts  for 
the  fact  that  the  orientation  of  the  missile  velocity  vector  is  not  necessarily 
along  the  instantaneous  LOS.  Since  missile  lateral  acceleration  is  devel¬ 
oped  normal  to  its  velocity  vector,  it  must  be  increased  by  the  indicated 

# 

factor  so  the  projection  normal  to  the  LOS  will  be  proportional  to  0. 

In  practice,  neither  0  or  0  are  directly  available  for  use  by 

It 

the  guidance  module.  If  the  guidance  law  is  operating  properly,  0  will  be 
small  for  all  but  the  last  fraction  of  a  second  before  intercept.  If  the 
variation  of  the  lead  angle  is  small,  then  the  factor  of  proportionality  can 
be  considered  constant  for  all  practical  purposes.  (This  is  also  consistent 
with  the  assumption  of  a  constant  closing  velocity  in  the  geometry  module). 
Some  missile  systems  use  the  gimbal  angle,  0h,  to  approximate  (0  -0) 
or  simply  enter  a  constant  which  is  representative  of  the  average  value.  ■ 
In  this  study,  a  constant  value  will  be  used  in  the  truth  model;  any  error 
in  this  value  can  be  accounted  for  as  an  effective  change  in  the  navigation 
ratio.  The  resulting  block  diagram  for  Guidance  Law  A  is  given  in  Fig. 

A.  6-1. 
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Figure  A.  6-1  Guidance  Law  A 


Kalman  Filter  for  Configurations  B,  C,  D  and  E  -  The  Kalman 
filter  is  implemented  in  the  non-rotating,  initial  LOS  coordinate  system. 
For  the  purpose  of  designating  the  filter,  it  is  assumed  that  both  the  target 
and  missile  maneuver  with  their  accelerations  normal  to  the  initial  LOS; 
the  resulting  plant  model  for  the  Kalman  filter  is  given  in  Fig.  A.  6-2. 


R- 8040a 
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i 

Figure  A.  6 -2  Kalman  Filter  Plant  Mode) 
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In  state -space  notation,  the  filter  model  equations  of  motion  are  given  by 
the  vector  differential  equation 


where 


x  (t)  =  Fx  (t)  +  u  (t)  +  b  (t) 


(A. 6-2) 
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Note  that  this  model  is  identical  to  the  system  truth  model  used  in  the 
target  and  geometry  modules  with  the  exception  that  target  and  missile 
acceleration  are  assumed  normal  to  the  orginal  LOS.  The  system  truth 
model  is  more  exact  in  that  the  missile  and  target  lateral  accelerations 
are  assumed  to  be  normal  to  their  respective  velocity  vectors,  which  in 
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general  do  not  lie  along  the  original  LOS.  Furthermore,  for  the  purpose 
of  mechanizing  the  filter,  it  is  assumed  that  the  autopilot  output  can 
be  measured  without  error. 

The  measurement  to  be  processed  by  the  Kalman  filter  is  the 
sampled  LOS  angle  6',  having  additive  independent  samples  of  noise  v^. 
The  latter  have  zero  mean  and  variance  given  by 

where  a  ,  a,  and  cr  represent  the  noise  components  defined  in  Section 
r  s 

A.  2.  In  terms  of  the  state  vector  x  in  Eq.  (A.  6-2),  the  LOS  angle 
measurement  is  expressed  as 


\  =h  (yxk+vk 


hT(t) 


L  tm 


W 


0  0 


] 


(A.  6-3) 


form 


The  discrete  Kalman  filter  mechanization  equations  have  the 


x  (t)  =  Fx  (t)  +  b  (t);  between  measurements  (A.  6-4) 


xk  =  x"  +  kk(zk  -  h1  (tfc)  x~ );  at  a  measurement  (A. 6-5) 


where  xk  denotes  the  solution  to  Eq.  (A.  6-4)  just  before  a  measurement 
is  processed.  The  gain  vector  kk  is  obtained  recursively  from  the 
matrix  covariance  equation  associated  with  the  Kalman  filter; 
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1 _ •> _  <- 


Mk  =  •Pk.1*T+S, 


pk  ■  \ 


■f 


eFW-r)Q[eF(At-r)]TdT 


PQ  =  E  jx(0)  x  (0)T}  $  =  eF^  (A. 6-6) 


The  gains  cannot  be  precomputed  and  stored  because  h  and  R  in  Eqs. 

(A.  6-6)  are  range  dependent;  this  requirement  is  responsible  for  most 
of  the  digital  computational  capability  required  by  the  guidance  system. 

The  mechanization  of  the  Kalman  filter  is  illustrated  in  Fig. 

^  ... . 

A.  6-3  where  and  at  are  the  respective  estimates  of  y^,  y^  and  a^.. 
Discrete  updates  of  the  estimates  occur  when  the  sampler  closes. 
Between  updates,  the  estimates  simply  propagate  according  to  the 
modelled  dynamics.  The  gains  k^,  k3k  and  k^  are  the  elements  ofk^ 
in  Eqs.  (A. 6-6). 

The  initial  covariance  matrix  P  (0)  used  in  the  computation  of 
the  Kalman  gains  contains  two  non -zero  elements.  The  initial  value  of 
its  third  diagonal  element,  p33,  is  the  variance  of  the  target  accelera¬ 
tion  at  the  initiation  of  the  terminal  phase.  It  will  be  assumed  that  the 
target  has  been  maneuvering  prior  to  the  start  of  terminal  phase  and 
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Figure  A.  6-3  Kalman  Filter  Mechanization  Diagram 


2 

the  initial  variance  is  equal  to  the  steady-state  maneuver  level,  j3 
(Eq.  A.  5-1).  The  initial  value  of  the  second  diagonal  element  of  P, 

Pgg,  is  equal  to  the  variance  of  the  relative  velocity  normal  to  the  initial 
LOS  at  the  initiation  of  the  terminal  phase.  From  Eq.  (A. 4-5)  of  Section 
A.  4,  it  can  be  shown  that  the  initial  value  of  y^  is  approximated  by 

h®  ~  -  <fce  vm  <A-6J7) 

for  small  values  of  initial  heading  error,  e. .  Therefore,  since  v  is 

n  m 

assumed  to  be  a  true  constant, 

P22(0)  =  “heV  <A'6-8> 

where  a^e  is  the  standard  deviation  of  initial  heading  error. 
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Controller  for  Guidance  LawsB,  C,  and  D-  The  control  input 
to  be  determined  is  the  commanded  missile  lateral  acceleration  in 
Fig.  A.  3»3  .  Continuous  control  will  be  assumed  in  deriving  the 
guidance  laws  and  it  is  desired  to  minimize  the  expected  square  of  the 
miss  distance  subject  to  a  penalty  function  on  the  total  control  energy. 
Therefore,  the  performance  index  to  be  minimized  is  given  by 

tf 

J  =  4  (tf)  +  y  J  aj,  (t)2  dt  (A.  6-9) 

*0 

subject  to  the  equations  of  motion  in  Eq.  (A.  6-2)  and  the 
autopilot  dynamics  modeled  by  the  first  order  transfer  function 


a 


!■ 


s  +  X 


m 


(A.  6-10) 


Note  that  the  saturation  in  Fig.  A.  3-3  is  neglected  in  Eq.  (A.  6-10),  as 
well  as  the  higher-order  autopilot  dynamics.  The  quantities  yd  (tf)  and 
y  are  respectively  the  terminal  miss  distance  at  intercept  time  t^  and 
the  weighting  on  control  effort.  The  solution  to  this  problem  is  called  an 
optimal  guidance  law.  By  invoking  the  separation  principle  (Ref.  10)  it 
is  known  that  the  control  will  be  of  the  form 

A  A 

ac=  clyd  +  ¥d  +  c3VcA  (A.6.11) 


The  indicated  control  gains,  which  minimize  the  performance  index  in 
Eq.  (A.  6-9), have  been  determined  by  Willems  in  Ref.  2  and  are 
repeated  below  for  convenience  (t  is  the  time  until  intercept): 


c 


n' 

go 


(A.  6-1 2) 
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n' 

c2=r 


go 


c3  =  n 


"M 


go 


7? 


go 


(A.  6-13) 


(A.  6-14) 


c4  = 


-n' 


.  x  t  .  i 

*m  go 
Am  go 


(A.  6-15) 


The  optical  navigation  ratio,  n',  is  given  by 


n'  = 


3y+  3  (i.;2^ 
K 


( x .  ;2  vgo) + .  2e- ^V) + ^  ^ .  y 


(A.  6-16) 

The  egression  for  n '  is  considerably  simplified  as  Xm  approaches  in¬ 
finity;  i.e., 


n‘ 


3t° 
go 

\n  =  °°  3y  +  tS 


(A.  6-17) 


go 


K  there  is  no  constraint  on  acceleration,  y  is  equal  to  zero  and  the 

resulting  navigation  ratio  from  Eq.  (A.  6-17)  is  constant  with  a  value  of 

3.  The  optimal  navigation  ratio  from  Eq.  (A.  6-16)  tends  to  increase  as 

t  approaches  zero.  If  y  is  equal  to  zero,  n '  will  approach  extremely 
go 

large  values  as  t  approaches  zero.  However,  for  a  non-zero  value  of 
go 

v ,  n'  will  first  increase  as  t  becomes  small;  after  a  point  it  will  begin 
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to  decrease  with  t  ,  and  finally  approach  zero  as  t  goes  to  zero. 

This  nonstationary  behavior  of  the  navigation  ratio  becomes  more  pro¬ 
nounced  as  X  becomes  small.  Increasing  n '  results  in  a  direct  scaling 
m 

of  all  four  control  gains,  which  in  turn  increases  the  autopilot  command, 
a  ,  to  compensate  for  the  dynamic  lag  represented  by  Eq.  (A.  6-10). 

V 


Given  this  set  of  optimal  linear  control  gains,  various  sub- 
optimal  approximations  can  be  made  to  simplify  the  computational  re¬ 
quirements.  If  v,  1/Xm,  cgand  04  are  arbitrarily  set  to  zero, classical 
proportional  navigation  with  n '  =  3  (Guidance  LawB)  is  the  resulting  con¬ 
trol  policy.  This  is  shown  by  differentiating  the  expression  for  LOS 
angle, 

yd 

6  *  =^-  (A.  6-18) 

mt 


to  obtain 


d  xmt  " : 

T~ 
mt 


d 


(A. 6-19) 


For  the  assumed  constant  closing  velocity  intercept, 


v„t 
c  go 


and  Eq.  (A.  6-19)  becomes 


6 


(A. 6-20) 


(A. 6-21) 
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Classical  proportional  navigation  is  defined  as 


a 


/  _ 


n  v 


(A.  6-22) 


Substituting  Eq.  (A.  6-21)  for  0  in  Eq.  (A.  6-22)  and  replacing  y^  and 
by  their  best  estimates,  obtained  from  the  Kalman  filter,  yields 


(A. 6-23) 


Comparing  Eq.  (A.  6-23)  with  Eqs.  (A.  6-11),  (A.  6-12),  (A.  6-13),  and 
(A.  6-17)  proves  the  equivalence  between  proportional  navigation  and  the 
optimal  guidance  law  under  the  conditions  1/Xm  =  y  -  c3  =  c^  =  0.  The 
navigation  ra*io  is  held  constant  in  Eq.  (A.  6-23).  Guidance  Law  C  is 
obtained  by  including  a  component  of  target  acceleration  in  the  formula¬ 
tion  of  the  autopilot  command  via  Cg.  Finally,  Guidance  Law  D  is 
realized  by  including  1  /\  4  0  in  the  computation  of  n ',  as  well  as  the 

term  involving  c4  in  Eq.  (A.  6-11).  In  all  cases  investigated  here  the 
weighting,  y ,  on  the  control  penalty  in  Eq.  (A. 6-9)  is  set  equal  to  zero. 

Optimal  Nonlinear  Controller  —  The  optimal  nonlinear  guidance 
law  is  based  u,,  n  the  same  system  model  used  for  the  optimal  linear  laws 
described  above;  however  its  performance  objective  is  different.  Name¬ 
ly,  an  acceleration  command  is  sought  such  that 

J  =  E  1 1 miss  distance  |  j  (A.  6-24) 

is  minimized,  subject  to  the  constraint 
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where  a _ is  the  maximum  airframe  acceleration.  Thus  the  limit  on 

max 

missile  acceleration  capability  is  incorporated  explicitly  into  the  prob¬ 
lem  formulation;  this  is  contrasted  with  the  quadratic  integral  penalty 
term  on  a'  (t)  in  Eq.  (A. 8-9). 

V 


The  solution  to  the  above  problem  is  derived  in  Refs.  7  and 
11.  The  resulting  optimal  nonlinear  guidance  law  has  the  same  struc¬ 
ture  as  the  linear  laws  discussed  previously;  however,  the  control  gains 
el’  c2*  c3  c4  are  comPu*ec*  differently.  The  latter  are  obtained  at 
each  control  computation  stage  as  a  result  of  the  following  steps: 


•  Determine  the  predicted  terminal  miss  distance 
based  on  the  Kalman  filter  estimate  of  the  state 
vector  at  the  current  control  computation  stage, 
and  the  dynamic  model  defined  in  Eqs.  (A.  6-2) 
and  (A.  6-10). 

•  Determine  the  value  of  commanded  acceleration, 
a^,  required  over  a  single  control  interval  to  null 
the  predicted  terminal  miss  distance,  neglecting 
the  constraint  in  Eq.  (A.  6-25). 

•  Determine  the  actual  acceleration  command,  ao, 
by  passing  a^  through  a  saturation  function 
that  satisfies  Eq.  (A.  6-25). 

To  carry  out  the  above  calculations,  we  first  combine  Eqs. 

(A.  6-2)  and  (A.  6-10)  into  the  form 


where 


x  (t)  =  Fx(t)  +  b  a„(t) 


"o  1  0  0  “ 

”  0 

0  0  1-1 

IV 

;b  = 

0 

0  0  -xt  o 

0 

0  0  0  -x_ 

A _ 

m 

m 

(A. 6-26) 


(A. 6-27) 


♦frhe  law  is  called  an  optimal  nonlinear  law  because  the  nonlinear  con¬ 
straint  on  ac  is  incorporated  in  the  problem  formulation. 


A-43 


73 


THE  ANALYTIC  SCIENCES  CORPORATION 

_ - _ -  <  *  _ />  . 


The  predicted  terminal  miss  based  on  a  state  estimate  at  time  t  is  given 

by 

yd(tf)  =  ST(tf,t)*(t)  (A.  8-28) 

where  S  (tf,t)  is  the  first  row  of  the  state  transition  matrix  corresponding 
to  P  in  Eq.  (A.  6-27).  The  effect  of  a  constant  command  ,  appliedoveran 
interval  At,  on  the  terminal  miss  distance  is  given  by 

/t+At 

*(t+At,  r)b  dT 
t 

^6c(t)a'  (A.  6-29) 

The  objective  is  to  select  a(.  such  that 

A  yd  (tf)  =  -  yd  (tf)  (A.  6-30) 

in  order  to  null  the  predicted  miss  distance.  Substituting  from  Eqs. 

(A.  6-28)  and  (A.  6-29)  leads  to  the  result 

ac  =  *  0T  (tf,t)x(t)- 

=  ”[c2^d  +  C2^d  +  c3  at  +  c4  ajg]  (A.  6-31) 


Carrying  out  the  computation  for  Cj  through  c4  we  find  that  they 
are  the  same  as  in  Eqs.  (A.  6-12)  through  (A. 6-15)  except  that  the  navi¬ 
gation  ratio  n'  is  given  by 
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4* 


ft- 


go 


At\ 

+  "T/ 


■f  \'V®' 

TTe 


\n 


(A.  6-32) 


The  values  of  n 'given  by  Eqs,  (A.  6-16)  and  (A. 6-32)  behave  quite  dif¬ 
ferently  as  functions  of  t  ,  as  demonstrated  in  Fig.  A.  6-4.  The  eaplana 
tion  for  this  is  that  the  nonlinear  law  attempts  to  completely  null  the 
predicted  miss  at  each  control  stage;  this  requires  much  larger  missile 
acceleration  early  in  the  trajectory  than  the  linear  laws,  which  tend  to 
reduce  the  predicted  miss  distance  more  gradually.  The  nonlinear  law 
has  a  smaller  gain  near  intercept  because  it  reduces  the  miss  distance 
faster  than  the  linear  law. 


Figure  A.  6-4  Comparison  of  Navigation  Ratios  for 
Optimal  Linear  and  Nonlinear  Laws 
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Autopilot  Command  -  The  autopilot  acceleration  command 
given  by  Eq.  (A.  6-11)  will  in  practice  be  contained  in  the  output  register 
of  a  digital  computer.  The  value  of  this  command  will  change  at  discrete 
intervals  which  are  nearly  coincident  with  the  LOS  measurements  pro¬ 
cessed  by  the  Kalman  filter.  In  practice  there  will  be  a  transport  lag 
between  the  LOS  measurement  and  the  update  of  the  control  signal  due  to 
the  computer  processing  time  requirements.  It  will  be  assumed  that  the 
time  delay  is  small  cc  pared  to  the  sampling  period  and  will  be  neg¬ 
lected.  Since  the  control  signal  changes  at  the  sample  instant  and  is 
constant  between  samples,  this  operation  is  equivalent  to  that  of  a  zero- 
order  hold.  Therefore,  the  autopilot  and  autopilot  model  are  driven  by 
a  stair-case  signal  which  is  the  sampled  and  held  acceleration  command 
given  by  Eq.  (A.  6-11). 


A.  7  COMPUTATIONAL  REQUIREMENTS 

The  computational  requirements  for  the  guidance  laws  dis¬ 
cussed  in  Section  A.  6  are  primarily  dictated  by  the  mechanization  equa¬ 
tions  for  the  Kalman  filter  in  Eq.  (A.  6-6).  Reference  15  provides  a 
guide  for  estimating  the  Kalman  filter  computational  requirements,  in 
terms  of  numbers  of  additions,  multiplications,  and  logical  operations 
required  to  complete  the  processing  of  each  measurement,  and  the 
associated  memory  capacity  needed.  Logical  operations  constitute  those 
involved  with  retrieving  variables  from  memory,  storing  in  memory, 
reading  instructions,  etc. ;  these  frequently  contribute  as  much  to  the  total 
computation  time  as  do  the  arithmetical  (add,  multiply,  etc.)  calculations. 

Using  the  egressions  derived  in  Ref.  15,  the  following  results 
are  obtained: 
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3 -State  Kalman  Filter  Computational  Requirements 

No.  of  Multiplications  a- 

No.  of  Additions  s- 

No.  of  Logical  Operation  Units  a* 

No.  of  Storage  Locations  a* 

Representative  execution  times  for  a  modem  minicomputer  are: 
multiplication  =  8  jusec;  addition  =  2  ^sec,  logical  operation  unit  =  l^sec. 
Applying  these  multiplication  factors  to  the  numbers  of  operations  item¬ 
ized  above,  and  allowing  for  the  use  of  two  3 -state  Kalman  filters  in  the 
3 -dimensional  engagement  situation,  the  total  execution  time  per 
measurement  cycle  and  storage  requirements  given  in  Table  A.  7-1  are 
obtained.  These  estimates  are  conservative  in  that  they  do  not  allow  for 
programming  efficiencies  that  can  be  gained  by  accounting  for  sparse 
matrices  (lots  of  zero  elements)  in  the  Kalman  filter  equations  (Eq. 

(A.  6-6)). 

TABLE  A.  7-1 

APPROXIMATE  KALMAN  FILTER 
COMPUTER  REQUIREMENTS: 

TWO  THREE -STATE  FILTERS 


125 

90 

2950 

Programming  - 

490 

Arrays 

60 

Operation 

Number  Required 

Computation  Time  (^sec) 

Addition 

180 

360 

Multiplication 

250 

Logical 

(Load,  Store,  etc.) 

5900' 

5900 
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APPENDIX  B 

THE  COVARIANCE  ANALYSIS  DESCRIBING 
FUNCTION  TECHNIQUE  (CADET) 


The  Covariance  Analysis  Describing  Function  Technique 
(CADET)  is  a  method  for  analytically  determining  the  statistical  pro¬ 
perties  of  a  nonlinear  system,  recently  developed  at  The  Analytic 
Sciences  Corporation  (Ref.  6).  The  principal  advantage  of  this  technique 
is  that  it  circumvents  monte  carlo  simulations,  thereby  achieving  sub¬ 
stantial  savings  in  computer  running  time.  We  first  motivate  the  dis¬ 
cussion  by  reviewing  covariance  analysis  methods  for  linear  systems; 
then  we  develop  an  analogous  procedure  (CADET)  for  the  nonlinear 
case. 


B.  1  COVARIANCE  ANALYSIS  FOR  LINEAR  SYSTEMS 

The  dynamics  of  a  linear  continuous  stochastic  system  can  be 
represented  by  a  first-order  vector -matrix  differential  equation  in  which 
x  (t)  is  the  system  state  yector  and  w  (t)  is  a  random  forcing  function, 

x (t)  =  F(t)  x(t)  +  w (t)  (B.l-1) 

Figure  B.l-1  illustrates  the  equation.  The  state  vector  is  composed  of 
any  set  of  quantities  sufficient  to  completely  describe  the  behavior  of  the 
system.  The  forcing  function  w(t)  represents  disturbances  as  well  as 
control  inputs,  that  may  act  upon  the  system.  It  can  be  assumed  without 
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R-230 


Figure  B.  1-1  Illustration  of  Continuous  Representation 
of  Linear  Dynamic  System  Equations 

loss  of  generality  that  the  system  states  and  the  forcing  function  are 
random  processes,  each  having  an  ensemble  average  value  of  zero;  i.e. , 
they  have  zero  mean.  In  what  follows  the  forcing  function  w(t)  is 
assumed  to  be  composed  of  elements  which  are  uncorrelated  in  time; 
that  is,  w(t)  is  ’’white*'  noise  having  the  spectral  density  matrix  Q(t) 
defined  by* 

E  [w(t)w(r)T]  =  Q(t)  6  (t-r) 

Under  the  above  conditions,  the  random  state  can  be  described 
in  terms  of  its  covariance  matrix  P  (t), 

P(t)  =  E  (x(t)x(t)T]  (B.l-2) 

The  equation  for  the  propagation  of  the  covariance  matrix  for  the  system 
described  by  Eq.  (B.l-1)  can  be  written  as  (Ref.  10). 

P(t)  =  F  (t)  P  (t)  +  P(t)F(t)T  +  Q(t)  (B.l-3) 


*E  denotes  ensemble  expectation ,  or  average  value;  T  denotes  matrix 
transpose. 
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The  diagonal  elements  of  F(t)  are  the  mean  square  values  of  the  state 
variables;  the  off-diagonal  elements  represent  the  amount  of  correlation 
between  different  state  variables. 

Equation  (B.l-8)  provides  a  direct  method  for  analyzing  the  sta¬ 
tistical  properties  of  x  (t).  This  is  to  be  contrasted  with  the  monte  carlo 
method  where  many  sample  trajectories  of  x(t)  are  calculated  from  com¬ 
puter-generated  random  noise,  or  random  numbers  in  the  case  of  a  digital 
computer.  If,  using  the  latter  technique,  m  such  trajectories  are  generat 
ed  using  Eq.  (B.  1*  1)— each  denoted  by  xfe(t),k  =1,  ...  ,  m— then  P  (t)  is 
given  approximately  by 

P(t)  2-  P(t)  =  4  £  5c(t)5c  (t)T  (B'1'4) 

k=l 

In  the  limit  as  m  approaches  infinity  we  have 

lim  P(t)  =  P(t)  (B.l-5) 

m  -♦  * 


Note  that  Eq.  (B.  1-3)  provides  an  exact  solution  for  P(t),  to 
within  computer  integration  accuracy,  whereas  the  monte  carlo  method 
yields  an  approximate  solution  for  a  finite  value  of  m.  Furthermore,  Eq. 
(B.  1-3)  need  be  solved  only  once  over  the  trajectory,  whereas  Eq. 
(B.l-1)  must  be  solved  many  times  using  the  monte  carlo  technique; 
consequently  the  direct  analytical  method  is  generally  the  most  efficient 
technique  for  analyzing  linear  systems.  Our  purpose  here  is  to  describe 
a  procedure  whereby  the  statistics  of  a  nonlinear  system  can  be  com¬ 


puted  approximately  using  a  recursion  relationship  similar  in  form  to 
Eq.  (B.l-3). 
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B.  2  COVARIANCE  ANALYSIS  FOR  NONLINEAR  SYSTEMS 

The  nonlinear  counterpart  of  Eq.  (B.l-1)  has  the  form 

x(t)  =  f(x(t))  +  w(t)  B.2-1) 

In  order  to  develop  a  covariance  analysis  method  similar  to  that  used 
for  liniar  systems  it  is  desirable  to  approximate ^(x)  in  Eq.  (B.2-1)  as 
a  linear  operation  onx(t).  In  particular,  we  assume  that  x(t)  is  com¬ 
posed  of  a  known  mean  x(t)  and  an  unbiased  additive  random  component  r(t), 

x(t)  =  x(t)  +  r(t)  (B.2-2) 

and  we  shall  seek  an  approximation  to  _f(x)  of  the  form 

y(t)  =  t(x(t) )  a  y(t)  +  Nf (t)  r(t)  (B. 2-3) 

where  Nf (t)  and  y(t)  are  to  be  specified.  The  structure  of  this  approxi¬ 
mation  is  illustrated  in  Fig.  B.2-1.  If  Eq.  (B. 2-3)  is  substituted  into 
Eq.  (B.2-1),  then  the  equations  of  motion  for  x(t)  become  linear; 

i(t)  =  £(t)  +  Nf(t)  (x(t)  -  x(t) )  +  w  (t)  (B.  2-4) 

Taking  the  expected  value  of  both  sites  of  Eq.  (B.  2-4)  yields 
an  equation  for  the  propagation  of  x(fc)~  namely, 

x(t)  =  i(t)  =  f(x(t))  (B.2-5) 

Subtracting  Eq.  (B. 2-5)  from  Eq.  (B.  2-4)  produces  an  equation  for 
r(t)  having  the  same  form  as  Eq,  (B.l-1); 

r(t)  =  Nf(t)r(t)  +  w(t)  (B.  2-6) 
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R-81R3 


Figure  B.  2-1  Block  Diagram  Interpretation  of  the 
Linearizing  Approximation 

Thus  by  analogy  with  Eq.  (B.  1-3), the  covariance  matrix  associated  with 
r(t)  —  also  referred  to  as  the  covariance  of  x— satisfies  the  differential 
equation 

P(t)  =  Nf(t)  m  +  P(t)Nf(t}T  +  Q(t)  (B.2-7) 

Equations  (B.2-5)  and  (B.  2-7)  together  provide  the  desired  analytical 
description  of  the  statistic s.of  x(t)— namely  its  mean  value  and  its 
covariance  matrix— assuming  that  y(t)  and  Nf(t)  are  known.  Subsequent 
paragraphs  discuss  the  approach  taken  in  CADET  for  determining  these 
quantities. 

Statistical  Linearization:  The  Scalar  Case --The  method  used 
in  CADET  for  approximating  f(x)  in  Eq.  (B.  2-1)  is  based  upon  statistical 
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linearization.  The  basic  principle  of  this  technique  is  conveniently 
illustrated  for  a  scalar  function,  f(x),  of  a  random  variable  x. 

Assume  that  f  (x)  is  to  be  approximated  by  the  linear  expres- 


f(x)  =  +  n^x 


(B.2-8) 


In  order  to  determine  appropriate  values  of  the  coefficients  n^  and  n^  we 
define  a  function  representation  error,  e,  of  the  form 


e  =  f (x)  -  nQ  -  nj  x 


(B.2-9) 


It  is  desirable  that  the  coefficients  be  chosen  so  that  e  is  small  in  some 
’•average”  sense;  the  procedure  used  to  accomplish  this  goal  is  to  mini¬ 
mize  the  mean  square  error,  E  [e2  ].  Thus,  forming 


E(e2)  =  E  [f2  +  n2  +  n2  x2  -  2nQf  -  2nxfx  +  2nonix] 


we  require  that 


(B.2-10) 


b  E  (e  ;  _  Q 
bn0 

bE  (e2)  _  n 


(B.2-11) 


which  are  necessary  conditions  for  a  minimum.  Solving  Eq.  (B.  2-11)  for 


no  “l  Produces 


Hq  =  f  -  n^x 


„  fx  -  fx 

1  =  ^  .2" 

X  -  X 


(B.2-12) 
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where  the  (A )  notation  denotes  expectation.  The  fact  that  these  values 
actually  minimize  the  mean  square  error  can  be  verified  by  examining 
the  second  derivative  of  E  (e  )  with  respect  to  n^  and  n^.  Substituting 
nQ  from  Eq.  (B.  2-12)  into  Eq.  (B. 2-8)  produces 

f(x)  =  f  +  n j  (x  -  x)  (B.2-13) 


Equation  (B.2-13)  constitutes  a  statistical  linearization  of  f 
about  the  mean  value  of  x  because  it  depends  upon  the  probability  density 
function  for  x,  p(x).  That  is, 

00 

f  =  J  f  (x)  p  (x)  dx 


J  f  (x)  x  p  (x)  dx 

-CO 


(B.  2-14) 


Observe  that  Eq.  (B.2-13)  is  in  the  desired  form  of  Eq.  (B. 2-1)  with  y 

A 

and  N|  identified  as  the  scalar  quantities,  f  and  n1  respectively. 

The  quantity  n^  in  Eq.  (B,  2-12)  is  usually  referred  to  as  the 
describing  function  gain.  Describing  functions  are  given  in  Ref.  12  for 
a  wide  variety  of  nonlinear  functions  with  gaussian  inputs.  From 

A 

knowledge  of  n^  and  f,  the  statistics  of  x  can  be  computed  using  Eqs. 
(B.2-5)  and  (B.  2-7)— hence  the  terminology  Covariance  Analysis 
Describing  Function  Technique  (CADET). 


A  more  common  method  of  linearizing  a  nonlinear  function  is 
the  use  of  a  first-order  Taylor  series  expansion 

f  (x)  =  f  (x)  +  -|L  I  ^  (x  -  x)  (B.  2-15) 

x=x 
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In  many  cases,  the  CADET  procedure  tends  to  be  the  more  accurate 
approach  from  a  statistical  point  of  view.  To  demonstrate  why  this  is 
so  consider  the  example  of  the  saturation  nonlinearity  in  Fig.  B.  2-2(a) 
and  assume  its  input  has  zero  mean.  If  f(x)  for  this  case  is  expanded  in 
a  Taylor  series  of  any  order  about  the  origin  (x  =  0),  we  obtain 

f  (x)  2;  x  (B.  2-16) 

The  effect  of  the  saturation  is  completely  lost  because  of  the  discontinuity 
in  the  first  derivative  of  f .  By  contrast,  if  statistical  linearization  is 
used,  we  have 


f(x)  2:  njX  (B.7-17) 

where  n^  is  the  describing  function  gain  defined  by 

r 

/  xf  (x)  p  (x)  dx 

nx  =  — ^ -  (B.2-18) 

j  x2  p  (x)  dx 


and  p  (x)  is  the  probability  density  function  for  x.  If  we  now  assume  that  x 
is  a  zero  mean  gaussian  random  variable,  then 

x2 

Pto  »  — e  2a  (B.2-19) 

\Z2ircr 

Substituting Eq.  (B.2-19)  into  Eq.  (B.2-18)  and  evaluating nj  for  the 
saturation  function  shown  in  Fig.  B.2-2(a),  we  obtain  (from  Ref.  12) 
the  result  shown  in  Fig.  B.  2-2(b).  It  can  be  seen  that  n^  is  a  function 


{a)  SATURATING  NONLINEARITY 


(b)  DESCRIBING  FUNCTION 


Figure  B.2-2 


The  Describing  Function  for  a  Saturation 
Nonlinearity  (Ref.  12) 


of  the  linear  part  of  f(x),  the  point  6  at  which  saturation  occurs,  and  the 
standard  deviation  of  x.  The  essential  feature  of  the  describing  function 


is  that  it  takes  into  account  the  probability  that  x  can  lie  within  the  sat 


uration  region. 


For  values  of  a  which  are  small  relative  to  fi  ,  the  probability 
of  saturation  is  low  and  n^  is  approximately  equal  to  1;  i.e. ,  fj  is 
approximately  equal  to  the  Ta$or  series  given  in  Eq.  (B.  2-16).  For 

larger  values  of  cr ,  n^  is  significantly  smaller  than  1  because  there  is  a 
higher  probability  of  saturation. 


As  a  result  of  the  above  discussion,  we  can  see  one  distinct 
advantp  that  statistical  linearization  has  over  the  Taylor  series  ex¬ 
pansion;  it  does  not  require  the  existence  of  derivatives  of  f(x).  Thus,  a 
large  number  of  nonlinearities— relays,  saturation,  threshold,  etc.— 
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can  be  treated  by  this  method  without  having  to  approximate  discon¬ 
tinuities  at  corners  in  f(x)  by  smooth  functions.  On  the  other  hand,  an 
apparent  disadvantage  of  the  method  is  that  the  probability  density  func- 

/V 

tion  for  x  must  be  known  in  order  to  compute  f  and  n^,  a  requirement 
that  does  not  exist  when  f  (x)  is  expanded  in  a  Taylor  series  about  its 
mean  value.  However,  as  discussed  in  Section  B.3,  approximations  can 
often  be  made  for  the  probability  density  function  such  that  the  resulting 
statistical  approximation  for  f  (x)  is  considerably  more  accurate  than  the 
Taylor  series,  from  a  statisiical  point  of  view. 


Statistical  Linearization— The  Vector  Case  -  We  now 
seek  a  linear  approximation  for  a  vector  function  f_(x)  of  a  vector 
random  variable  x,  having  probability  density  function  p  (x) .  Following 
the  statistical  approximation  technique  outlined  for  the  scalar  case,  we 
propose  to  approximate  _f(x)  by  the  egression 


_f(x)  a*  a  +  Nf  x 


(B .  2-20) 


where  a  and  Nf,  are  a  vector  and  a  matrix  to  be  determined.  Defining 
the  error 


e  =  l(x)  -  a  -  Nf  x 


(B.  2-21) 


a  and  are  chosen  so  that  the  quantity 

J  =  E  [eT  Ae  ] 


(B.  2-22) 


is  minimized  for  some  symmetric  positive  semidefinite  matrix  A.  Sub¬ 
stituting  Eq.  (B.  2-21)  into  Eq.  (B.  2-22)  and  setting  the  partial  derivative 
of  J  with  respect  to  the  elements  of  a  equal  to  zero,  we  obtain 


j 
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Therefore,  a  is  given  by 


a  =f.  (x)  -  Nf  x 


(B.2-24) 


Substituting  a  from  Eq.  (B.2-24)  into  J  and  taking  the  partial  deriva¬ 
tive  with  respect  to  the  elements  of  yields 

E  ^A[N^xxT  +  (f(x)  -!(x))xT]]  =0  (B.2-25) 

whe'.e 


~  /V 

X  =  X  -  X 


Solving  Eq.  (B.2-25)  produces 


/\ 


T  T  _-l 


Nf  =  [Is  -I  X  ]P 


(B.2-26) 


where  P  is  the  covariance  matrix  of  x.  Observe  that  both  a  and  Nf  as 
given  by  Eqs.  (B.  2-24)  and  (B.2-26)  are  independent  of  the  weighting 
matrix  A;  hence,  they  provide  a  generalized  minimum  mean  square 
error  approximation  to  f. 


Upon  substitution  from  Eq.  (B.2-24)  into  Eq.  (B.2-20)  we 

obtain 

f  (x)  =  ?  (x)  +  Nj  (x  -  x)  (B.  2-27) 


which  is  the  form  specified  in  Eq.  (B.  2-3).  The  quantity  Nf  is  called 
the  describing  function  matrix;  as  expected,  Nf  reduces  to  n<  in  Eq. 

(B.  2-12)  when  xandj[  are  scalars.  Because  is  potentially  a  function 
of  both  the  mean  and  covariance  of  x,  denoted  by  Nf  (P,  x),  equations 
(B.  2-5)  and  (B,  2-7)  become  a  set  of  coupled  nonlinear  differential 
equations; 
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|(t)  -£(t) 

P(t)  =  Nf  (P,  x)  P  +  PNf  (P,  x)T  +  Q(t)  (B.2-28) 


B.3  COMPUTATION  OF  DESCRIBING  FUNCTIONS 

In  order  to  carry  out  the  integration  of  Eq.  (B.  2-28),  it  is 
necessary  to  compute  y(t)  and  N^(t).  From  the  discussion  in  Section 
B.2,  we  know  that 

7(t)  =f(t) 


Nf(t)  =  [fxT  -  hT]pJ 


(B.3-1) 


where  x  and  P  are  the  mean  and  covariances  of  x(t),  respectively  and 


i  =  J  f(x)p(x)dx 


/N 

T 

fx  = 


/  j[  (x)  xT  p(x)  dx 


(B.3-2) 


Thus  x  and  P  can  be  continuously  evaluated  if  the  probability  density 
function  p(x)  in  Eq.  (B.3-2)  Is  known. 

Generally  speaking,  it  is  not  practical  to  analytically  evaluate 
the  probability  density  function  for  the  state  of  a  nonlinear  dynamical 
system;  therefore  an  approximate  form  of  p  (x)  must  be  obtained.  One 
frequently  used  assumption  is  that  x  is  gaussion: 
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p(x)M2»Det(P)rn/2exp  [-  |(x-i)'IV1  (*  -  x>  1  (B.3-3) 


This  is  based  upon  the  fact  that  the  states  of  the  nonlinear  system  are 
the  result  of  integrating  the  quantity  f:  (x)  f  w.  Therefore  x(t)  is  a  linear 
superposition  of  past  values  of  random  quantities,  which  will  tend  to  be 
gaussion  even  through  f  (x)  and  w  may  not  be  nongaussian.  The  latter 
assertion  is  based  on  qualitative  application  of  the  central  limit  theorum 
(Ref.  8).  If  Eqs.  (B. 3-2)  and  (B.3-3)  are  combined  with  Eq.  (B.2-28), 

A 

x(t)  and  P  (t)  can  be  evaluated  numerically  on  a  digital  computer. 

From  the  above  discussion,  it  is  clear  that  more  computation 
is  required  in  applying  covariance  analysis  techniques  to  nonlinear  sys¬ 
tems  than  for  linear  systems,  because  of  the  procedure  used  to  evaluate 
f  and  Nf.  Consequently,  efficient  means  for  calculating  the  right  hand 
sides  of  Eq.  (B.  3-2)  are  desired.  One  useful  simplification  arises  when 
the  nonlinear  system  consists  of  a  few  single-state-input  single-output 
nonlinearities,  such  as  that  illustrated  in  Fig.  B.3-1.  In  this  case, 
making  the  gaussion  assumption  described  above,  each  nonlinear  opera¬ 
tion  is  replaced  by  its  mean  value  and  its  describing  function  gain  as 
indicated  in  Fig.  B.  2-1,  both  of  which  are  scalars  computed  according 
to  Eq.  (B.  2-12).  Then  the  first-order  vector  and  matrix  differential 
equations  in  Eqs.  (B.  2-28)  can  be  formulated  by  inspection  of  the  result¬ 
ing  linear  system  block  diagram.  Thus  the  CADET  equations  can  be 
constructed  using  existing  tabulated  describing  functions  for  single-input 
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where  the  term  Fx  describes  the  dynamics  of  the  linear  part  of  the  sys¬ 
tem.  Applying  Eq.  (B.3-4),  we  have 


(B.3-6) 


The  term  bf^/bXj  is  simply  the  describing  function  gain  for  f.,  as 
determined  by  the  statistics  of  its  input,  x^;  thus  Nf  is  determined  by 
the  dynamics  (F)  of  the  linear  system  and  the  single-input  describing 
gain.  The  extension  of  the  above  argument  to  the  case  when  several 
single -state -input  nonlinearities  exist  is  straightforward. 


B.4  MIXED  CONTINUOUS -DISCRETE  SYSTEMS 

Preceeding  sections  of  this  appendix  have  treated  continuous 
nonlinear  systems;  i.e.,  those  governed  by  differential  equations. 
However,  in  many  practical  applications,  the  system  may  include  a 
digital  computer  whose  input  and  output  are  expressed  in  terms  of 
difference  equations,  as  illustrated  in  Fig.  B.4-1.  Such  a  structure 
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arises  in  missile  guidance  systems  where  digital  control  laws  are  used 
to  command  the  missile's  maneuvering  acceleration.  In  this  section, 
equations  are  developed  for  propagating  the  mean  and  covariance  of  a 
nonlinear,  mixed  continuous-discrete  system. 

The  equations  of  motion  for  a  system  of  the  type  shown  in 
Fig.  B.4-1  are  written  in  mixed  differential-difference  equation  format. 
First  of  all,  between  sampling  intervals  the  digital  computer  is  effectively 
idle  and  the  continuous  part  of  the  system  satisfies  an  equation  of  the 
form 

x,(t)  =  f(x,(t))+ £(^(0)  + w(t)  \ 

^(t)  .  o  j  ‘k£t<W  !B-4‘1) 

where  x  (t)  refers  to  the  continuously  varying  states  in  the  system,  and 
x^(t)  is  a  collection  of  digital  states  (e.g.,  states  in  the  digital  computer) 
which  remain  unchanged  during  the  sampling  interval.  Assuming  that 
the  composite  state  vector,  x,  defined  by 


RANDOM  DISCRETE  R-II2J 

INPUT  0YNAM1CS 


Figure  B.4-1  An  Example  of  a  Mixed  Continuous 
Discrete  System 
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is  gaussian  with  its  mean  and  covariance  matrix  known  at  time  t^, 
then  its  statistics  during  the  sampling  interval  can  be  calculated  using 


d 


dt  xd(t) 


=  [|(xc(t))+g(2d)  ];  tk  st  <t( 


Nf(t)  Ng' 

P  (t)  =  P(t)  +  P(t) 

0  0 


Nf(t)T  0 


T 

Ng  0 


Oft)  o 


0  0 


tkst<tk+l 


(B.  4-2) 


where  and  are  the  describing  function  gain  matrixes  associated 
with  f^  and  £  respectively.  Observe  that  N  remains  constant  throughout 
the  sampling  interval  because  x^  is  itself  constant. 

Now,  at  a  sampling  interval,  tne  digital  computer  performs  a 
calculation  which  can  be  represented  as  a  difference  equation  of  the  form 


^dvvk+l  ' 


■^  k+1  '»  ^d'  k+1 


-k+1 


(B.4-3) 


where  the  superscript;  (-)  denotes  the  solution  to  Eq.  (B.4-1)  just  before 
the  sampling  instant  and  (+)  denotes  the  new  values  of  the  state  variables 
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just  after  a  sampling  instant.  The  vector  represents  a  random 
quantity  that  can  enter  the  digital  calculation  as  a  result  of  instrument 
measurement  noise.  It  is  assumed  that  w^+j  has  zero  mean  and 
covariance  matrix  Q^.  Observe  that  in  Eq.  (B.4-3)  x,  remains  un¬ 
changed  because  only  the  digital  states  can  change  instantaneously  in 
time. 


Because  the  mean  and  covariance  of  x  and  x.  at  t,  ..  are 

■— c  — d  k+1 

known  from  Eq.  (B.4-2),  the  describing  function  gain  matrix  cor¬ 
responding  to  h  in  Eq.  (B„  4-3)  can  be  evaluated.  Thus  we  can  rewrite 
Eq.  (B.4-3)  approximately  as 

2c  (tk+l+)  =  2 c  <W> 


2d  <tk+l+)  ~  ^  +  “h  (5(tk+r>  -  i(tk+1"))  + 

(B.4-41 

From  Eq.  (B.4-4)  it  follows  that  the  mean  and  covariance  of  the  system 
states  just  after  the  sample -time  are  given  by 


2A+1+> 

2d<‘k+l+> 


2<W> 

^2c(tk+x'>-  2A+r>)_ 


f 

H* 

>_  + 

II 

_I !  0 
_l  _ 

p'‘k+r> 

"i ! 

I  T 

--'K 

+ 

0  1  0 

1 

1 

~l  ”  ™ 

1 

L 

0  I 

(B.4-5) 

After  evaluating  Eqs. (B.4-5),  x(t^+j+)  and  P  (tk+j+)  are  the  initial 
conditions  for  propagating  the  mean  and  covariance  over  the  next 
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sampling  interval  using Eqs .  (B .  4-2).  Thus  by  alternately  evaluating  the 
differential  and  difference  propagation  equations,  the  mean  and  covariance 
matrix  of  the  mixed  continuous -discrete  nonlinear  system  can  be  evaluated. 


B.  5  DESCRIBING  FUNCTIONS  FOR  SATURATION  AND  SINUSOIDAL 
NONLINE  ARTHES 


Describing  functions  for  the  saturation  and  sinusoidal  non- 
liniarities  in  Fig.  2.1-2  are  calculated  in  Ref.  12  for  gaussian  inputs. 
These  are  shown  graphically  in  Figs.  B. 5-1  and B. 5-2.  The  exact 
analytical  expressions  are  as  follows: 


Saturation 


n(<r)  =  m  2PI  (£-)  -1 


PI 


Sinusoidal 


n  (or)  =  Mme 


-m2(r2/2 


i 


s 


(B.5-1) 


(B.5-2) 


Equation  (B.5-1)  contains  a  term,  PI  -j  ,  which  is  the  integral  of  the 

gaussian  density  function,  sometimes  called  a  probability  integral. 

The  latter  cannot  be  expressed  in  closed  form.  However,  an  approxi- 

-5 

mate  expression  which  is  accurate  to  one  part  in  10  is  given  by  (Ref.  14) 
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b  =  0.33267 
at  =  0.4361836 


a„  =  -0.1201676 
2  (B.5-3) 

a„  =  0.9372980 


Equations  (B.5-1)  through  (B.5-3)  are  used  to  calculate  the  describing 
function  gains  in  the  CADET  simulation. 


R-5156a 


Figure  B.5-1 


Describing  Function  for  Saturation 
with  a  Gaussian  Input 
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Figure  B.5-2  Describing  Function  for  a  Sinusoid 
with  a  Gaussion  Input 
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